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INTRODUCTION 


The  original  purpose  of  this  research  was  to  investigate  three  outstanding  problems 
associated  with  the  use  of  ocean  bottom  seismographs  and  hydrophones: 

1.  mechanical  coupling  between  an  OBS  and  the  ocean  floor  for  both  horizontal  and 
vertical  motion; 

2.  hydrodynamic  coupling  between  an  OBS  and  the  near  bottom  water,  both  at 
seismic  amplitudes  and  frequencies  and  at  those  associated  with  near  bottom 
currents; 

3.  the  relative  advantages  and  disadvantages  of  seismometers  versus  hydrophones, 
singly  or  in  arrays,  for  various  hydroacoustic  and  seismic  observations. 

Because  of  a  reduced  funding  level  the  third  problem  was  not  studied  in  detail. 

The  problem  of  coupling  an  OBS  to  soft,  ocean  bottom  sediments  has  been  of  con¬ 
cern  for  some  time.  Recognizing  this  problem,  ONR  sponsored  the  Lopez  Island  OBS 
intercomparison  experiment  (Sutton  et  al,  1980  and  1981a).  This  experiment  produced  a 
great  deal  of  quantitative  information  and  provided  a  considerable  amount  of  insight  on 
directions  to  go  for  improving  design.  Taking  advantage  of  knowledge  gained  at  Lopez 
Island,  two  follow-up  OBS  coupling  tests  were  conducted  by  personnel  of  Hawaii  Insti¬ 
tute  of  Geophysics  in  cooperation  with  RAI  in  Keehi  Lagoon,  Oahu  in  (Duennebier  et  al, 
1981,  Sutton  et  al,  1981b,  Tuthill  et  al,  1981).  In  these  tests,  six  different  HIG  opera¬ 
tional  and  developmental  OBS  configurations  were  compared  using  calibrated  mechani¬ 
cal  transients,  sine  wave  and  white  noise  drivers,  and  explosions.  A  large  amount  of 
excellent  data  was  generated  during  the  Keehi  Lagoon  tests.  We  helped  complete  the 
analyses  of  this  data  and  prepare  the  results  for  publication,  emphasizing  the  quantita¬ 
tive  aspects  of  the  mechanical  coupling  tests  and  the  amount  of  distortion  expected  for 
different  OBS  configurations. 

The  least  understood,  or  studied,  aspect  of  OBS  coupling  is  the  effect  of  the  near 
bottom  water,  both  in  production  of  noise  by  interaction  of  currents  and  in  signal  distor¬ 
tion  through  seismic/hydroacoustic  forces  applied  directly  to  the  body  of  an  OBS. 
Definite  differences  of  "opinion"  regarding  the  importance  of  water  coupling  became  evi¬ 
dent  during  and  subsequent  to  the  Lopez  Island  experiment.  We  have  developed  an 
appropriate  theory  and  compared  it  with  OBS  data  available  at  HIG  in  order  to  obtain  a 
reliable  estimate  of  the  importance  of  water  coupling  as  a  function  of  OBS  configuration. 
Laboratory  tests  employing  a  mass  suspended  in  water  on  a  shaker  table  confirmed 
theoretical  predictions. 

The  paper  "Optimum  Design  of  Ocean  Bottom  Seismometers"  by  G.  H.  Sutton  and 
F.  K.  Duennebier,  which  is  submitted  to  Marine  Geophysical  Researches,  summarizes  the 
results  of  this  research.  An  invited  paper  "Ocean  Bottom  Seismology:  History  and 
Current  Status"  by  G.  H.  Sutton,  presented  at  the  5ACLANT  Ocean  Seismo-Acoustics 
Symposium  in  June  1985  and  published  in  the  Proceeding  volume  (1986,  T.  Akal  and  J. 
M.  Berkson,  eda.)  includes  some  of  these  results  and  a  summary  of  the  current  status  of 
OBS  technology.  Both  of  these  papers  are  included  in  this  Final  Report.  "A  Comment 
on  ’the  Recovery  of  True  Particle  Motion  from  Three-Componet  Ocean  Bottom 
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’  Seismometer  Data’  by  J.  D.  Garmany"  by  G.  H.  Sutton,  F.  K.  Duennebier,  and  G.  J. 

,  Fryer  is  also  included  in  this  report.  The  shaker  table  used  in  the  laboratory  tests  was 

described  in  a  paper  by  F.  K.  Duennebier,  G.  H.  Sutton,  D.  Harris,  and  D.  Byrne,  "A 
k  Simple  Shaker  Table  for  Seismometer  Calibration",  Marine  Geophysical  Research,  6,  P. 

y  311-328,  1984. 
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Abstract 

Ocean  bottom  seismometers  (OBS)  have  been  widely  used  during  the  past  decade  to  collect 
seismic  data  for  determination  of  the  structure  of  the  oceanic  lithosphere,  stress  patterns  in 
regions  of  earthquake  activity,  and  geo-acoustic  parameters  of  the  ocean  floor.  Data  quality 
from  these  experiments  has  often  been  disappointing  because  of  poor  signal  quality  and  high 
noise  levels.  Many  of  these  problems  result  from  motion  of  the  OBS  package  that  is  decoupled 
horn  motion  of  the  ocean  floor.  These  coupling  problems  are  more  serious  in  the  ocean  than 
on  land  because  of  the  low  shear  strengths  of  most  ocean  sediments.  In  this  paper  we  continue 
to  develop  the  theory  of  coupling  of  OBSs  to  soft  sediments  and  arrive  at  results  suggesting 
that  OBS  packages  should  be  designed  with:  (l)  the  minimum  mass  possible,  (2)  radius  of  area 
in  contact  with  the  sediment  proportional  to  the  cube  root  of  the  mass,  and  the  maximum 
radius  less  than  1/4  of  the  shear  wavelength,  (3)  density  of  the  OBS  approximately  that  of  the 
sediment,  (4)  a  low  profile  and  a  small  vertical  cross  section  with  the  water,  and  (5)  low  density 
gradients,  and  maximum  symmetry  about  the  vertical  axis.  Agreement  of  the  theory  with  test 
data  is  good;  most  deviations  are  reasonable,  given  limitations  of  the  theory  and  experiments. 
The  theory  also  suggests  that  the  coupling  frequency,  the  frequency  above  which  the  OBS  does 
not  follow  the  motion  of  the  sediment,  is  directly  proportional  to  the  sediment  shear  velocity. 


1  INTRODUCTION 

Seismometers  are  instruments  designed  to  measure  the  motion  of  the  material  they  are  attached 
to.  A  well-designed  seismometer  will  detect  this  motion  with  high  fidelity  and  low  noise  over 
a  broad  range  of  amplitudes  and  frequencies,  modifying  the  motion  only  by  a  known  transfer 
function.  These  characteristics  are  relatively  easy  to  achieve  on  land  by  careful  placement  of 
seismometers  in  low  noise  environments  on  materials  with  high  shear  strengths.  [Krohn  (1984)] 


1 


I 


discusses  problems  associated  with  coupling  of  exploration  geophones  to  different  soils  for  fre¬ 
quencies  above  100  Hz.  There  are  few  locations  on  the  ocean  bottom  with  high  shear  strengths, 
and  finding  a  low  noise  site  is  a  matter  of  chance.  High-quality  data  can  be  obtained  from  the 
ocean  bottom  by:  (1)  avoiding  experiments  in  soft  sediments,  (2)  using  only  hydrophones,  which 
are  much  less  sensitive  to  coupling  problems  than  seismometers  but  can  not  directly  detect  shear 
arrivals  [Koelsch  and  Purdy  (1979)j,  (3)  avoiding  the  low-shear  sediments  entirely  by  emplacing 
seismometers  in  boreholes,  which  is  very  expensive  [Duennebier  et.al.  (1981), Byrne  et.aL  (1983)], 
(4)  trying  to  determine  the  change  in  the  transfer  function  caused  by  coupling  problems 
[Zelikoritz  and  Prothero  (1981)],  or  (5)  designing  OBSs  that  will  not  be  sensitive  to  coupling  prob¬ 
lems.  Although  several  attempts  have  been  made  to  design  OBSs  to  minimize  coupling  prob¬ 
lems  [Duschenes  et.al.  (1981),Trehu  and  Solomon  (1981), Byrne  et.al.  (1983)],  the  attempts  have 
suffered  from  the  lack  of  a  coupling  theory  that  defines  the  parameters  that  are  critical  to  OBS 
design. 

In  this  paper  we  present  theoretical  arguments  that  define  the  parameters  important  to  design 
of  OBSs.  Results  are  compared  with  data  obtained  from  the  Lopez  Island  OBS  Intercomparison 
Experiment  [Sutton  et.al.  (1981b)],  two  experiments  conducted  in  Keehi  Lagoon,  Oahu,  and  labo¬ 
ratory  experiments  on  a  shaker  table.  Two  theories  are  presented  to  explain  the  motions  of  OBSs 
on  soft  sediments,  one  (the  spring-mass  theory)  is  appropriate  for  low  frequencies  and  sediments 
that  are  treated  as  an  elastic  solid,  but  it  is  difficult  to  extrapolate  this  theory  to  high  frequencies. 
The  other  theory  (the  fluid  theory)  is  appropriate  at  liigh  frequencies  and  when  the  sediment  can 
be  treated  as  a  fluid.  These  theories  agree  well  with  each  other  except  for  one  term,  the  amount  of 
sediment  entrained  in  the  motion  of  the  OBS. 

The  spring-mass  theory  was  originally  developed  by  structural  engineers  to  analyze  the  vibra¬ 
tions  of  foundations.  It  has  been  used  by  seismologists  to  model  the  behavior  of  seismometers  in 
swamps,  and  was  applied  to  the  OBS  case  by  including  buoyancy  [Sutton  et.al.  (1981a)].  In  the 
previous  paper,  we  had  no  theoretical  basis  to  determine  how  much  water  and  sediment  are  in- 
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eluded  in  the  motion  of  the  OBS  when  it  moves  relative  to  the  sediment.  In  this  paper,  the  theory 
is  developed  to  include  the  entrained  water  and  sediment.  The  fluid  theory  is  a  result  of  studies  of 
the  motions  of  bodies  immersed  in  fluids  (Bachelor  (1967)].  The  motions  of  spheres  and  disks  in 
fluids  is  studied  to  obtain  equations  for  the  motions  of  OBSs  in  fluids.  The  resulting  formulations 
are  compared  with  experimental  field  data  for  coupling  of  vertical  motion,  and  with  laboratory 
data  for  both  vertical  and  horizontal  motion.  Each  theory  is  valid  when  used  to  predict  the  motion 
of  OBSs  in  different  situations.  The  theories  agree  well  with  the  data,  and  clearly  imply  guidelines 
for  design  of  OBSs  for  optimum  coupling  to  the  bottom,  despite  problems  of  not  understanding 
the  motions  in  all  cases.  In  addition  to  the  design  constraints  implied  by  the  theory,  other  design 
constraints,  such  as  small  cross  section  in  the  water  and  symmetry  about  the  vertical  axis,  are  also 
discussed  to  show  ways  of  minimizing  cross  coupling  and  coupling  of  water  motion  to  motion  of  the 
OBS. 

2  THEORY 

2.1  Spring-Mass  Theory 

A  seismograph  sitting  on  a  compliant  ocean  floor  has  six  degrees  of  freedom  of  motion:  translation 
and  rotation  in  the  vertical  and  two  orthogonal  directions.  For  our  purposes,  we  shall  ignore  rotation 
about  the  vertical  axis,  and  concentrate  on  translation  in  the  vertical  direction.  Each  degree  of 
freedom  can  support  a  mode  of  oscillation  with  a  characteristic  frequency  and  damping.  In  general, 
the  modes  will  be  coupled.  For  a  well-designed  instrument,  however,  only  horizontal  translation  and 
rocking  about  the  orthogonad  horizontal  axis  may  be  strongly  coupled.  Radial  symmetry  about  the 
vertical  axis,  which  does  not  exist,  for  example,  with  the  commonly  used  tripod  support,  is  required 
for  horizontal  coupled  modes  to  be  independent  of  azimuth.  The  vertical  modes  are  decoupled  with 
a  good  tripod  mount. 

Forces  impinge  on  the  OBS  from  both  the  sediment  and  the  water.  Differential  horizontal 
motion  across  the  water-sediment  boundary,  both  from  seismic  waves  and  from  bottom  currents, 
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can  produce  distortion,  cross  coupling,  and  noise. 

Much  of  the  theory  required  for  the  analysis  of  dynamic  coupling  between  an  OBS  and  the 
bottom  has  been  developed  by  structural  engineers  for  the  design  of  foundations  resting  on  soils. 
[Hsieh  ( 1962))  provided  a  clear  summary  of  this  work  along  with  references  to  important  ear¬ 
lier  work.  [Safar  (1978)]  applied  the  theory  to  investigate  the  minimization  of  distortion  caused 
by  coupling  between  the  ground  and  vertical  geophones.  (He  uses  the  dynamical  interaction 
between  a  number  of  geophones  to  improve  the  transient  response  to  seismic  signals.)  These 
studies  do  not  include  the  effects  of  the  water.  [Sutton  et.al.  (1981a)],  [Sutton  et.aL  (1981b)], 
(Zelikovitz  and  Prothero  (1981)],  and  [Trehu  (1985a)]  consider  the  coupling  and  seismic  response 
of  an  OBS  surrounded  by  water.  This  paper  is  an  extension  of  that  work. 

First,  we  review  the  equations  appropriate  for  the  motion  of  an  instrument  on  a  uniform  elastic 
half-space  in  the  absence  of  water,  and  then  include  the  effect  of  the  water.  Only  vertical  motion 
will  be  considered  in  detail.  Horizontal  coupling  and  the  effect  of  layering  or  gradients  within  the 
bottom  will  be  considered  briefly. 

The  detailed  analysis  of  the  dynamic  response  of  a  rigid  body  resting  on  a  solid  half-space  is 
quite  complex.  In  general,  the  frequency  response  depends  upon  the  instrument  mass,  and  the 
size,  shape,  and  stress  distribution  of  the  contact  area,  and  the  shear  velocity,  density,  dissipation, 
and  Poisson’s  ratio  versus  depth  in  the  half-space.  [Hsieh  (1962)]  provides  equations  appropriate 
for  a  rigid  circular  disk  on  a  homogeneous,  isotropic  elastic  medium  for  vertical  and  horizontal 
translation  and  for  rotation  about  vertical  and  horizontal  axes.  The  coefficients  of  the  solutions 
depend  upon  Poisson’s  ratio  and  the  ratio  of  the  contact  radius  to  the  wavelength  of  shear  waves 
in  the  medium.  These  coefficients  can  be  modified  for  different  shapes  and  distributions  of  contact 
pressure.  (Safar  (1978)]  summarized  equivalent  results  for  vertical  motion  only.  [Luco  (1974)] 
calculated  impedance  functions  for  rigid  foundations  on  a  layered  medium  and  concluded  that,  in 
contrast  with  a  uniform  half-space,  a  layered  medium  shows  stronger  frequency  dependence  and 
reduced  radiation  damping  (higher  Q  resonances)  at  low  frequencies. 
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Typical  ocean  bottom  sediment  can  have  a  shear  velocity,  v,,  less  than  50  m/s  and,  perhaps, 
as  low  as  10  m/s;  compressional  velocity,  vp,  is  near  that  for  water,  1.5  km/s;  and  density,  p,, 
ranges  between  about  1300  and  1300  kg/m3.  Thus,  for  the  relations  used  in  this  paper,  we  shall 
assume  Poisson’s  ratio,  a  =  0.5  (nearly  fluid)  and  p  —  1570  kg/m3.  The  equations  presented  are 
appropriate  only  for  shear  wavelengths,  A„  greater  than  4  times  the  bearing  radius,  r.  This  can  be 
a  serious  restriction  since,  for  example,  at  v ,  =  20  m/s  and  20  Hz,  the  radius  should  be  less  than 
0.25  in. 

The  response  of  the  seismograph  frame,  /,  to  a  vertical  seismic  input  signal,  2,  in  the  absence 
of  water  and  assuming  a  negligible  mass  of  sediment  is  displaced  by  the  instrument  is  given  by 
[Safar  (1978)]  as 


I/Z  =  (KS/Q)  +  u>c2]/[SJ  +  (uicS/Q)  +  u,'1 

]  (1) 

where  =  [K(a)/ M,\1^  =  coupling  frequency 

(2) 

Q  =  [ K(a)Mi]l/2/R(a )  =  quality  factor 

(3) 

Mi  —  OBS  mass 

and  5  =  the  Laplace  transform  variable. 

The  dynamical  spring  constant,  K,  and  the  damping  coefficient,  R,  are 

both  functions  of  the 

dimensionless  frequency  given  by, 

a  =  wr/v,  =  2xr/A, 

and,  therefore,  wc  and  Q  in  equation  1  are  also  functions  of  frequency. 

For  a  circular  disk  on  an  elastic  half  space, 

K  -  vlp,rFx(a)  and  R  =  v,p,r3Fj(a)/a. 

(4) 

For  <r  =  0.5  and  0  <  a  <  1.5, 

F\  2;  8.0  -  2.0 a2  and  Fj  ~  6.9a  [Hsieh  (1962)]. 

(5) 
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Fi  anil  h\  have  different  values  for  different  shaped  OBSs,  and  for  horizontal  translation  or  rotation 
about  horizontal  or  vertical  axes. 

In  Figure  1,  the  approximate  values  for  F\  and  F2  given  by  equations  5  are  shown  to  agree  well  FIG  I 
with  values  obtained  from  detailed  calculations  by  Bycroft  [Bycroft  (1956), Arnold  et.al.  (1955)]. 
Substituting  equations  4  and  5  in  equations  2  and  3  and  solving  for  u  =  aic  we  obtain 


=  [K0/M'}1/2  (6) 

and 

Q  =  1/2  h  =  acb/6.9  =  0.41b/(6  +  2)1/2  (7) 

where  K0  -  8p,v?r  -  static  spring  constant,  <r  ~  0.5  [Hsieh  (1962)]  (8) 

A/'  =  M,  +  2 p,r3  —  effective  mass,  instrument  plus  sediment  (9) 

ac  =  wcr/u,  =  dimensionless  resonant  frequency  (10) 

b  -  M,/p,r3  =  dimensionless  shape  factor  (11) 


Some  algebra  also  shows  that  a*  =  8/(b  +  2)  and,  since  equation  5  is  valid  only  for  values  of  a 
less  than  about  1.5.  this  requires  that  5  >  1.56,  i.e.,  the  instrument  mass,  M,  >  1.56p,r3.  Figure 
2  shows  the  increase  of  effective  mass,  Ai\  and  fraction  of  critical  damping,  hc,  with  increase  in  FIG  2 
dimensionless  resonant  frequency,  ac. 

Because  of  the  dependence  of  Kq  and  M'  on  r,  for  given  (disk- shaped)  instrument  mass  and 
sediment  parameters,  the  coupling  frequency  exhibits  a  maximum  value  with  increasing  r.  This 
can  be  seen  in  Figure  3,  which  is  a  plot  of  coupling  frequency  versus  bearing  radius  for  different  FIG  3 
values  of  instrument  mass,  Mi,  from  Equation  6,  8,  and  9.  Assuming  that  the  maximum  coupling 
frequency  is  desirable,  the  optimum  radius  and  maximum  coupling  frequency,  also  obtained  from 
Equations  6,  8,  and  9  is: 

^  =  (Mi/4p.)l,i  (12) 

fcmax  =  0.292v.(p,/Mi)l/3.  (13) 
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At  the  optimum  coupling  condition:  6  =  4,  ac  =  1.16  and  Q  =  0.67.  The  increase  in  ropt  and 
decrease  in  maximum  coupling  frequency,  fc  mox,  with  increasing  instrument  mass  obtained  from 
Equations  12  and  13,  are  shown  in  Figure  4. 

2.2  Fluid  Theory 

The  above  analysis  does  not  include  any  effects  of  the  water  surrounding  the  OBS  above  the 
sediment.  Before  including  this  effect,  and  its  complicated  boundary  conditions,  we  will  study  the 
motion  of  an  OBS  surrounded  only  by  a  fluid. 

The  water  affects  the  seismic  response  of  the  instrument  in  two  ways:  1)  a  buoyancy  effect  that 
reduces  the  effective  force  of  the  seismic  signal  on  the  instrument  and  depends  upon  the  volume 
of  water  (and  sediment)  displaced;  and,  2)  an  added  virtual  mass  that  depends  upon  the  cross- 
sectional  area  and  shape  of  the  instrument.  [Bachelor  (1967)]  gives  the  equation  for  the  motion  of 
a  rigid  sphere,  u,>  in  response  to  the  motion  of  an  acoustic  wave  in  fluid,  vwi  for  wavelengths  large 
compared  to  the  diameter  of  the  sphere: 

vi/vw  =  3/(2(M</Mw)  +  1)  (14) 

where  Mj  and  Mw  are  the  masses  of  the  sphere  and  displaced  water,  respectively.  Figure  5  is  a 
plot  of  Equation  14  for  Mi/Mw  between  0  and  3.  This  indicates  amplification  approaching  a  factor 
of  3  for  luw  sphere  densities  and  attenuation  for  high  densities.  When  the  densities  of  the  sphere 
and  water  are  equal  (neutral  buoyancy)  the  motions  are  equal.  Many  operational  OBSs  have  a 
large  vertical  gradient  in  density,  with  a  dense  lower  support  structure  and  buoyant  superstructure. 
Based  on  the  above  analysis,  horizontal  signals  are  expected  to  produce  strong  torques  in  such 
instruments. 

The  virtual  mass  Af„,  for  a  sphere  immersed  in  water  is  Mw/ 2;  that  for  a  cylinder  is  Mw;  and 
that  for  a  disk  is  (8/3)p„,rJ  [Bachelor  (1967)].  The  motion  of  a  disk  in  water,  when  the  relative 
motion  is  normal  to  the  plane  of  the  disk,  is  given  by: 


FIG 


FIG 


.  1  +  8r/3rt 

Vi/Vw  —  7  77  “ 

Pi/ Pw  +  8r/3ict 


(15) 


where  pi  and  pw  are  the  disk  and  water  densities,  t  is  the  disk  thickness,  and  r  is  the  disk  radius. 
Equation  15  is  plotted  in  Figure  6  showing  v,/vw  as  a  function  of  r/t  for  several  values  of  pi/pw.  FIG  6 
This  motion  is  much  different  than  that  of  a  sphere  in  that  the  radius  of  the  sphere  has  no  effect 
on  its  motion,  but  the  ’’flatness”  of  the  disk  is  a  very  important  parameter  in  its  motion.  This 
change  in  response  with  shape  is  related  to  the  ease  with  which  water  can  move  past  the  body  in 
the  direction  of  motion.  The  flatter  the  disk,  the  more  difficult  it  is  for  water  to  move  past  the 
disk,  and  the  larger  the  mass  of  water  (virtual  mass)  that  moves  with  the  disk.  In  the  OBS  case, 
one  wants  the  OBS  to  move  with  the  sediment,  not  with  the  water,  thus  a  large  virtual  sediment 
mass  (relative  to  Mi)  and  small  virtual  water  mass  are  desired. 

On  a  fluid-fluid  boundary,  we  can  approximate  the  motion  of  a  disk  between  two  fluids  by  simply 
dividing  the  virtual  mass  into  two  equal  parts  with  the  appropriate  densities:  M„  =  4 (p,  +  p^,)r3  fZ. 

We  have  no  rigorous  basis  for  doing  this,  and  this  formula  does  not  account  for  energy  radiated 
away  along  the  boundary.  We  expect,  however,  that  this  value  of  Af„  is  a  good  approximation  for 
small  motions.  With  this  definition  for  virtual  mass  of  a  disk,  the  velocity  of  the  disk  relative  to  a 
fluid-fluid  interface  for  motion  perpendicular  to  the  interface  is  : 

i  _  g&a  +  4r/3vt 

This  equation  is  the  same  as  Equation  15  if  pw  =  p,- 

Here  we  note  that  if  pw  =  0,  we  should  have  the  same  result  as  for  Equation  1  when  K o  =  0, 
or  in  the  high  frequency  limit  (<v  2>  uic).  This  is  not  the  case  because  Equation  1  approaches  zero 
when  u>  >  uic,  and  Equation  16  is  a  constant  given  by  ( +- 1)-1  when  pw  =  0.  When  the  spring 
constant  is  set  to  zero  in  Equation  l,  I jZ  —  [(MiS/R)  +  l]-1,  which  is  a  low  pass  filter  function. 

While  both  equations  are  valid  in  their  range  of  interest,  they  yield  very  different  results  in  these 
cases.  Looking  more  closely  at  the  differences,  the  fluid  theory  includes  no  radiation  of  energy  away 
from  the  body,  while  the  spring-mass  theory  does.  The  increase  in  virtual  mass  generally  observed 
near  a  boundary  where  flow  is  restricted  is  also  not  included  [Milne-Thompson  (1950)].  We  expect 
that  Equation  1  yields  the  correct  results  on  an  interface  when  K o  0,  and  we  now  include  the 
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effects  of  water  above  the  interface  in  the  equations. 

2.3  Spring-Mass  theory  with  water 

Since  the  bottom  sediments  are  ’’nearly”  fluid,  in  this  paper  we  shall  use  the  fluid  dynamic  relations 
as  a  first  approximation  to  the  virtual  mass  of  the  water  (in  addition  to  the  term,  already  included 
in  Equation  9  for  the  sediment  side  of  the  boundary  obtained  from  the  dynamical  spring  constant). 
Note  that  the  virtual  mass  is  a  function  of  shape,  and  it  can  change  with  the  direction  of  motion. 

The  vertical  component  of  the  equation  of  motion  for  the  instrument  frame  including  the  effect 
of  water  can  be  written  in  Laplace  transform  notation  as 

(Mi  +  M„)S2I  +  RS(I  -Z)  +  K(J  -  Z)  =  (M„  +  M.  +  M„)S2Z  (17) 


where  Mw  and  M,  represent  the  mass  of  water  and  sediment,  respectively,  displaced  by  the  instru¬ 
ment,  and  M„  is  its  virtual  mass  of  water  entrained  in  the  motion.  Note  that  for  a  sphere  in  the 
absence  of  sediment,  i.e.  R  =  K  =  M,  =  0  and  Af„  =  Mw/2,  Equation  17  reduces  to  Equation 
14  that  gives  the  response  of  a  sphere  to  an  acoustic  wave  in  a  fluid.  From  Equation  17  we  can 
obtain  an  expression  for  the  response  of  the  OBS  frame  to  a  vertical  seismic  input,  in  the  presence 
of  water, 


a(S2  +  vc'S/aQ'  +  vc'2fa) 

1  (5’ WWW2) 


(18) 


where 


a  =  (Mw  +  M,  +  Af„)/Af* 


(19) 


and  we’  and  Q1  are  shape-dependent  parameters,  given  by  Equations  6  to  11  for  disk-shaped  OBS 
footpads,  with  the  instrument  mass,  M,,  replaced  by  the  dynamic  mass: 


M*  =  (M{  +  Mv). 

Equation  18  is  similar  to  that  given  by  (Sutton  et.al.  (1981a)],  with  their  coupling  coefficient,  C, 
replaced  by  (l  -  a).  However,  the  dynamic  mass  of  sediment  included  in  ‘he  earlier  definition  of  C 
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is  not  represented  in  a.  Sutton  et  al.  defined  (1  -  C),  equivalent  to  our  a ,  as 


t 


I 


(1-C)  = 


Mw  +  Ma  +  M„  +  A// 
Ml  +  m; 


(20) 


where  Ml  represents  the  sediment  that  moves  with  the  OBS.  Comparing  Equations  10  and  20,  we 
see  that  a  <  (1  -  C).  Although  an  effective  sediment  mass  (2 p.r3,  arising  from  the  dynamical 
spring  constant  K,  for  a  disk-shaped  area  in  contact  with  the  sediment)  is  included  in  the  resonant 
frequency  and  damping  terms  through  Equations  6  and  7,  apparently  it  is  not  required  as  an 
additional  virtual  mass  to  that  of  the  water,  M„.  In  the  absence  of  water,  Mw  =  Mv  =  0,  and  with 
negligible  displaced  sediment,  M,  =  0,  a  becomes  zero  and  Equation  18  reduces  to  Equation  1.  We 
note,  however,  that  this  equation  does  not  yield  the  desired  result  (Equation  15)  as  the  interface 
disappears  when  K  approaches  zero,  and  p,  approaches  p„  for  a  disk-shaped  OBS. 

[Sutton  et.al.  (1981a)J  added  a  buoyancy  term  to  the  restoring  force  supplied  by  the  elasticity 
of  the  sediment: 

Kt  =  K  +  (p,  -  pm)*r2g 


where  g  is  the  acceleration  of  gravity.  It  is  unlikely  that  the  latter  term  will  be  important  in  practical 
situations,  as  shown  below.  The  ratio  of  the  first  term  (A")  to  the  second  is,  using  Equation  8  and 
taking  Ko  -  K, 

K  _  8 v)p,T 

( P ,  ~  Pv,)*r2g  I0trr2(p,  -  pw) ' 

For  a  sediment  with  v,  -  10 m/s  and  p,  =  1300  kg/m 3  this  ratio  has  a  value  of  110/r,  thus  Kj-  is 
not  affected  by  buoyancy  unless  r  is  very  large. 

From  Equation  18  we  see  that  the  response  of  the  OBS  frame  to  a  seismic  input  signal  can 
be  defined  with  three  parameters:  w';  Q';  and  a.  These,  in  turn,  depend  upon  the  shape  and 
mass  of  the  instrument  package  and  the  elastic  properties  of  the  bottom  sediment.  The  transfer 
function  is  1  for  all  frequencies  when  the  coupling  coefficient,  a,  equals  1  i.e.,  neutral  buoyancy 
(Mi  =  Mw  +  M,  in  Equation  19)  and  no  net  dynamic  force.  By  substituting  ju/  for  5,  the  Laplace 
transform  variable,  in  Equation  18,  we  can  separate  the  real  and  imaginary  parts,  and  obtain  the 
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Note  that  the  above  equations  and  18,19,21,  and  22  do  not  depend  explicitly  on  OBS  shape. 

The  curves  in  Figures  (7)  and  (8)  are  for  the  approximation  that  the  coefficients  Q'  (=  1/2/t') 
and  w'  are  constant.  When  response  curves  using  constant  coefficients  are  compared  with  those 
obtained  by  using  coefficients  that  vary  with  frequency  according  to  Equations  4  and  5,  the  response 
is  within  1  dB  amplitude  and  5"  phase  difference  for  <*»  <  l.5w'.  For  this  reason,  we  ignore  the 
frequency  variation  of  u and  Q'  in  the  remainder  of  this  paper. 

As  indicated  above,  the  parameters  u Q',  and  a  may  be  calculated  using  equations  6  to  11 
with  Mi  replaced  by  M*  (for  disk-shaped  area  in  contact  with  the  sediment)  and  Equation  19  for  a 
particular  instrument  design  and  sediment,  as  Af„  is  strongly  dependent  on  the  shape  of  the  OBS 
in  the  water.  These  parameters  can  also  be  determined  experimentally  from  an  in  situ  calibration, 
as  discussed  below. 


2.4  Transient  tests 


A  force  can  be  applied  to  an  OBS  by  the  release  of  a  float  [Sutton,  et,al.,  1981a],  or,  alternatively, 
by  use  of  an  internal  shaker  [Zelikovitz  and  Prothero  (1981)].  Such  calibration  transients  have  been 


applied  to  both  operational  and  prototype  instruments.  The  motion  of  the  OBS  caused  by  a  force, 
T ,  on  the  OBS  is  given  by  the  following  equation: 


[M*(S2+wc'S/Q'WJ)l  V  ' 

Equation  23  can  be  obtained  from  equation  17  by  replacing  the  right  hand  side  by  T  and  setting 
2  =  0.  Using  Equation  23  or  its  Laplace  transform,  the  dynamic  mass  M’  and  the  coupling 
frequency,  w',  and  Q'  can  be  obtained  by  matching  the  observed  spectra  or  the  time  transients 
(after  correcting  for  the  response  of  the  geophone)  to  theoretical  response  functions  ">r  their  time 
equivalents,  for  a  known  calibration  force,  F.  The  coupling  coefficient,  a,  can  be  calculated  using 
Equation  19.  The  use  of  Equation  23  should  be  valid  independent  of  the  exact  nature  of  virtual 
mass,  Mv,  and  whether  or  not  the  sediment  contributes  to  that  term.  The  amplitude  and  phase 
functions  obtained  from  Equation  23  for  a  float  transient  test  (F  =  F/S)  are  given  by: 


(fo  =  fan 


-i  lw? 


2/i'uj'ui 


3  Experimental  Data 

3.1  Laboratory  test 


A  Laboratory  test  was  conducted  using  the  shaker  table  facility  at  the  Hawaii  Institue  of 
GeophysicsfDuennebier  et.al.  (1984)]  to  verify  the  responses  of  a  sphere  to  seismic  waves  predicted 
by  equations  21  and  22.  In  this  test  a  platform  isolated  from  ground  vibrations  by  a  weak  spring  is 
excited  with  motion  having  a  near-white  spectrum  in  a  vertical  or  horizontal  direction.  The  motion 
is  detected  by  a  geophone  rigidly  mounted  to  the  platform,  and  by  an  identical  geophone  mounted 
in  a  sphere  suspended  from  the  platform  by  a  spring  representing  the  sediment  (Figure  9).  The  FIG  9 
transfer  function  of  the  suspended  geophone  (T  in  the  figure)  relative  to  the  reference  geophone 
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(R)  yields  the  coupling  response,  and  should  be  predicted  by  Equations  21,  and  22  within  the 
experimental  errors. 

Two  tests  are  discussed  here,  one  for  vertical  motion  using  a  buoyant  sphere,  and  one  for 
horizontal  motion  using  a  heavy  sphere,  both  motions  were  measured  with  the  spheres  in  cur,  and 
immersed  in  water.  Data  on  the  spheres  and  test  results  are  given  in  Table  1.  The  low  density  T  1 
sphere  was  used  in  the  vertical  test  to  avoid  problems  of  non-linearity  in  the  suspension.  By  using 
a  density  of  about  O.bgm/cm3,  the  spring  extension  in  air,  is  the  same  as  it  is  in  water,  although 
the  spring  is  mounted  on  top  of  the  sphere  in  air  and  on  the  bottom  in  water.  The  amplitude  and 
phase  results  from  the  tests  in  air,  digitized  from  the  output  of  a  spectrum  analyser,  are  shown  in 
Figure  10  together  with  suitable  solutions  of  equations  21  and  22  with  a  =  0.  Note  that  the  theory  FIG  10 
and  data  are  in  excellent  agreement. 

When  the  spheres  are  immersed  in  water,  two  effects  are  expected:  1)  the  frequency  of  resonance 
should  decrease  because  of  the  added  virtual  mass;  and  2)  at  high  frequencies  the  amplitude  response 
should  become  constant  and  the  phase  response  should  return  to  zero  because  of  the  buoyancy  effect. 

As  seen  in  Figure  11,  each  of  these  effects  is  observed.  For  the  vertical  case,  the  change  in  a  is  FIG  11 
as  predicted  using  the  constants  computed  from  Equation  19  and  Table  1  within  the  experimental 
errors  (the  spectra  were  digitized  from  3”  by  4”  photographs  taken  from  the  Spectrum  Analyser). 

The  changes  observed  in  the  horizontal  test  are  as  predicted  except  that  the  value  of  a  predicted 
is  about  twice  as  large  as  that  observed  from  the  data,  as  shown  by  the  frequency  of  the  spectral 
low  at  about  11  Hz,  and  the  lower  than  expected  response  at  high  frequencies.  The  reason  for  this 
discrepancy  is  not  known,  and  is  puzzling  since  the  change  in  coupling  frequency  from  air  to  water  is 
very  close  to  that  predicted  (8.3  Hz  to  7.05  Hz).  It  is  likely  that  this  problem  results  from  the  finite 
size  of  the  water  tank  and  shallow  water  depth  used  to  immerse  the  sphere.  This  failure  is  serious 
in  that  proper  estimates  of  the  value  of  a  are  vital  in  order  to  correct  transfer  functions  obtained 
by  internal  calibration  devices  to  the  sesimic  transfer  function.  It  appears  that  a  is  sensitive  to 
the  geometry  of  the  instrument  and  to  the  relationship  to  nearby  interfaces.  If  a  is  not  accurately 
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■ 

Vertical  Test 

Horizontal  Test 

Sphere  diameter,  cm 

12.7 

10.2 

Sphere  volume,  cm3 

r'“  1070 

549 

Sphere  mass,  gm 

454 

850 

Sphere  density,  gm/cm3 

0.424 

1.55 

Mass  displaced  water,  gm 

1070 

549 

Dynamic  mass,  [M,  +  Af, ),  gm 

989 

1125 

Coupling  frequency,  Hi,  observed  in  air 

7.2 

8.6 

Coupling  frequency,  Hi,  observed  in  water 

4.7 

7.2 

Calculated  a  in  water,  from  Eqn.  19 

1.62 

0.73 

Observed  a  in  water,  from  Figure  1 1 

1.61 

0.37 

Observed  h,  from  Figures  10,  11 

0.03 

0.03 

Table  1:  TABLE  1.  HIG  Shaker  Table  Test  Parameters 


estimated  then  large  errors  in  predicted  amplitude  and  phase  response  can  result.  Uncertainties  at 
frequencies  above  the  coupling  frequency  resulting  from  the  variation  of  the  coefficients  w'  and  Q' 
are  also  likely. 


3.2  Field  Tests 


Several  major  field  experiments  have  been  conducted  within  the  past  few  years  to  test  the  coupling 
characteristics  of  existing  OBSs.  The  first  (and  largest)  was  the  Lopez  Island  OBS  Intercomparison 
Test  conducted  during  the  summer  of  1978,  during  which  the  characteristics  of  twelve  OBSs  were 
studied  [Sutton  et.al.  (1981b)|.  A  second  large  cooperative  test,  conducted  in  Brest,  Prance,  was 
mainly  for  European  groups  [Snoek  et.al.(1982)].  Three  smaller  tests  were  conducted  by  the  Hawaii 
Institute  of  Geophysics  in  Keehi  Lagoon,  Oahu,  in  1981,  1982  and  1985,  to  test  modifications  made 
to  instruments  suggested  by  the  results  of  Lopez.  Trehu  (1985a)  reported  result  of  tests  on  the 
OBS  used  by  the  U.S.  Geological  Survey.  The  primary  test  in  each  case  was  the  float  transient  test 
explained  earlier.  In  addition,  a  series  of  shots  were  fired  to  supply  seismir  ources  in  some  of  the 
tests.  In  this  paper  we  restrict  the  discussion  to  the  transient  test  results. 

Instrument  parameters  and  results  from  the  transient  tests  are  given  in  Table  II.  T  2 
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The  instrument  masses,  equivalent  radii  of  the  area  in  contact  with  the  bottom,  and  mass  of 
the  water  displaced  by  the  OBS  are  ail  quantities  dependent  only  on  instrument  design.  Qth, 
ropt>  fcmax i  Af„,  M*,  and  a  are  obtained  from  instrument  parameters  and  an  assumed  sediment 
density  of  p,  =  1570&J /m3  and  sediment  velocity  of  20  m/s.  fc,  and  v,  are  obtained  from 
the  response  of  the  vertical  geophone  in  each  OBS  to  a  vertical  float  transient  test,  and  equations 
presented  above.  An  explanation  of  each  parameter  in  the  table  is  given  below: 

NAME:  Name  of  OBS  tested  using  the  acronyms  from  the  LOPEZ  Experiment.  Names  ending  in 
”L”  are  from  the  LOPEZ  Experiment,  names  ending  in  ”1”  and  ”2”  are  from  tests  in  Keehi 
Lagoon,  Oahu,  in  1081  and  1982.  Instruments  are  listed  in  groups  by  experiment,  and  by 
increasing  dynamic  mass  within  experiments. 

A/,-:  Instrument  mass  in  kg. 

A/*,:  Mass  of  water  displaced  by  the  instrument. 

r:  Equivalent  radius  of  the  area  in  contact  with  the  sediment. 

fc:  Coupling  resonance  frequency  (and  Q)  are  measured  by  comparison  of  spectra  from  transient 
tests  with  theoretical  spectra  from  equation  24.  An  example  is  shown  in  Figure  12. 

Qot,,:  Quality  factor  of  the  coupling  frequency  resonance  peak  measured  from  spectral  analysis  of 
the  transient  test  data. 

Qth'-  Quality  factor  of  the  coupling  frequency  resonance  peak  from  Equation  7. 
v,:  Shear  velocity  of  the  sediments  calculated  from  Equation  8. 

r«pt:  Optimum  radius  for  area  in  contact  with  the  sediment  calculated  from  Equation  12  with  A/; 
replaced  by  M * . 

Jc  max'  Maximum  coupling  frequency  possible  (at  r  =  r^t)  from  Equation  13. 


FIG  12 
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M„  Virtual  mass  of  water  entrained  in  the  motion  of  the  OBS.  This  mass  depends  on  the  shape 
of  the  OBS,  and  is  estimated  by  summing  the  contributions  of  different  parts  of  the  OBS, 
including  the  water  above  the  area  in  contact  with  the  sediment  (4pwr3/3),  the  water  around 
any  spheres  (Mw/2),  and  the  water  around  any  (horizontal)  cylinders  (Mw). 

M’:  The  dynamic  mass  of  the  OBS  when  emersed  in  water  (M*  =  M,  +  Mv). 

a:  The  coupling  parameter  defined  by  Equation  19. 

From  Equations  6,  8,  and  9  we  see  that  the  quantity  fc/v„  the  coupling  frequency  divided  by 
the  sediment  shear  velocity,  should  be  a  simple  function  of  the  sediment  density,  the  OBS  bearing 
radius,  and  Af*.  As  normal  sediment  densities  occur  only  within  bounds  from  about  1300  to  1700 
kg/m3,  this  function  is  mainly  determined  by  the  dynamic  mass  (M*  =  +  Af„)  and  the  bearing 

radius  (r)  of  the  OBS.  Thus,  for  a  given  OBS,  a  measurement  of  fc  should  be  directly  proportional 
to  the  shear  velocity.  A  plot  of  fc/v,  calculated  from  Equation  6  (using  an  appropriate  dynamic 
mass  and  radius  from  Table  II)  versus  fc  measured  from  transient  test  spectral  data  (Table  II)  is 
shown  in  Figure  13.  Each  point  represents  the  average  of  several  coupling  frequency  measurements  FIG  13 
at  either  Lopez  Island  or  Keehi  Lagoon.  The  shear  velocity  of  the  sediment  beneath  each  OBS  can 
be  estimated  for  each  point  by  dividing  the  observed  fc  by  the  predicted  fc/v,.  Although  there  is 
considerable  scatter,  note  that  all  but  a  few  of  the  27  points  lie  between  15  and  30  m/s,  and,  as  two 
different  sites  are  represented  and  each  OBS  was  in  a  slightly  different  location,  much  of  the  scatter 
is  probably  the  result  of  real  variations  in  shear  velocity.  There  is  no  obvious  correlation  of  shear 
velocity  with  instrument  mass  or  bearing  radius,  suggesting  that  the  theory  successfully  accounts 
for  changes  in  these  parameters.  Note  that  (with  the  exception  of  the  Plate  Standard,  having  a 
mass  of  140  kg  when  tested  at  Lopez  Island)  all  OBSs  with  observed  coupling  frequencies  above 
15  Hz  have  masses  less  than  50  kg,  and  all  OBSs  with  coupling  frequencies  less  than  10  Hz  have 
masses  greater  than  100  kg.  OBSs  with  masses  that  differ  by  more  than  an  order  of  magnitude  yield 
estimates  of  shear  velocity  that  agree  within  a  factor  of  2,  suggesting  some  validity  for  the  theory. 
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The  same  data  are  presented  in  Figure  14  in  a  different  way.  In  Figure  14,  using  the  relationships 
in  Table  II,  the  ratio  of  observed  fc  to  fc  max  is  plotted  versus  the  ratio  of  actual  radius  to  optimum 
radius.  The  theoretical  curve  is  shown  for  comparison.  As  for  the  preceeding  figure,  much  of  the 
scatter  probably  results  from  real  variations  in  sediment  shear  velocity. 

Q'  of  the  coupling  resonance  at  /e  as  predicted  by  Equation  7  is  plotted  versus  observed  Q'  in 
Figure  15.  It  is  noted  that  the  predicted  Q'  is  lower  them  the  observed  Q'  in  most  cases,  which 
means  that  most  OBSs  "bounce”  at  the  coupling  frequency  more  than  is  predicted  from  the  half¬ 
space  theory.  Note  that  all  but  one  observed  Q‘  value  above  2.  are  from  large-mass  OBSs,  and  all 
but  one  Q'  value  less  than  2  are  from  small-mass  OBSs.  Thus,  although  the  theory  does  not  match 
the  observed  values  of  Q',  the  desirability  of  a  low  Q'  coupling  resonance  still  leads  us  to  the  same 
conclusion:  that  small  mass  OBSs  perform  better  than  large  mass  OBSs.  The  reason  the  theory 
does  so  poorly  in  modeling  Q'  is  not  known;  however,  there  are  several  possibilities.  The  theory 
assumes  that  the  sediment  is  a  homogeneous  half  space,  but  velocity  generally  increases  with  depth. 
This  gradient  should  be  "felt”  by  the  more  massive  OBSs  with  the  greater  bearing  radii  and  with 
the  lower  coupling  frequencies.  Energy  would  not  radiate  away  from  the  system  as  fast  as  predicted 
by  the  half-space  theory,  in  agreement  with  the  layered  theory  of  [Luco  ( 1974)].  The  smaller  OBSs 
would  not  feel  as  deeply. 


4  OBS  DESIGN  PARAMETERS 

On  the  basis  of  the  theory  and  observations  discussed  above,  an  ocean  bottom  seismometer  »hou . 
be  designed  with  the  following  constraints: 

1.  minimum  mass, 


FIG  H 


FIG  15 
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2.  radius  in  contact  with  the  ocean  floor  given  by  equation  12  (with  M,  replaced  by  A/*),  and 
radius  less  than  1/4  of  the  shear  wavelengths,  and 

3.  density  approximately  that  of  the  sediment. 

In  addition,  the  following  constraints  will  minimize  rocking  modes,  cross  coupling,  and  noise 
coupled  to  the  OBS  by  motion  of  the  water: 

4.  low  height- to-base  area  ratio,  low  moment  of  inertia  about  horizontal  axes,  and  small  ( smooth ) 
cross  section  with  the  water,  and 

5.  low  density  gradients,  and  maximum  symmetry  about  the  vertical  axis. 

Each  of  the  above  constraints  is  discussed  below: 

1).  As  data  near  and  above  the  coupling  frequency  are  potentially  distorted  by  motion  of  the 
OBS  relative  to  the  sediment,  it  is  advisable  to  keep  this  frequency  as  high  as  possible  by  minimizing 
the  OBS  mass.  On  land  (and  in  the  ocean  where  there  is  no  sediment  cover)  this  is  generally  not 
a  problem  because  shear  velocities  are  usually  high  enough  to  move  the  coupling  frequency  above 
the  range  of  interest. 

Note  that  while  this  analysis  is  done  for  vertical  coupling,  a  similar  analysis  can  be  done  for 
horizontal  coupling.  As  there  are  more  degrees  of  freedom  in  horizontal  coupling,  the  situation 
is  more  complicated.  Horizontal  transient  test  data  suggest  that  the  coupling  problem  is  more 
severe  in  the  horizontal  direction  than  in  the  vertical,  and  that  values  of  Q'  a r»  higher  and  coupling 
frequencies  are  lower.  The  small-mass  (moment  of  inertia)  guideline  still  holds,  however,  since  the 
same  parameters  are  involved.  Note  also  that  Q‘  tends  to  increase  with  increasing  mass  (Figure 
15),  and  that  low  Q’  at  coupling  resonances  is  also  desirable. 

If  mass  is  not  minimized,  the  response  of  the  OBS  in  the  seismic  frequency  band  will  potentially 
vary  with  shear  velocity  of  the  near-bottom  sediments.  Although,  as  mentioned  earlier,  it  is  possible 
to  make  in  situ  calibrations  using  internally  generated  transients,  the  relationship  between  the 
seismic  response  and  transient  test  response  will  be  in  doubt  (for  instance,  if  the  OBS  is  not 
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vertical).  According  to  the  theory,  if  the  OBS  has  the  same  density  as  the  surrounding  material, 
it  should  not  be  affected  by  coupling  as  long  as  the  seismic  wavelengths  are  much  larger  than  the 
OBS  (see  3  below). 

2) .  If  the  equivalent  radius  of  the  area  in  contact  with  the  sediment  is  not  optimal,  then  the 
coupling  frequency  will  be  lower  than  the  best  possible  value  for  that  mass  of  OBS  (iW'J  and 
sediment  density  (Figure  3).  OBSs  with  relatively  small  radii  with  respect  to  their  mass,  (such  as 
I/DGOJL  and  UWL  in  Table  II)  have  low  coupling  frequencies,  and  PS1  and  PS2,  with  relatively 
large  radii,  also  show  low  coupling  frequencies.  Only  three  OBS  tests  (PSL,  OSU12,  and  OSU22) 
show  coupling  frequencies  significantly  higher  than  expected  from  the  OBS  characteristics  and  those 
of  the  sediments.  The  reason  is  unknown,  but  they  may  have  been  resting  on  stilfer  sediment  or 
rocks.  Trehu  (1985a)  reports  results  of  in-situ  transient  tests  on  OBSs  of  differing  mass  and  base 
radius.  These  tests  are  for  configurations  with  ratios  of  mass-to-base  radius  greater  than  those 
covered  by  the  second  order  approximation  to  K  and  R  given  in  equations  4  and  5.  Her  observed 
resonant  frequencies  are  coiupatable  with  a  zero  order  theory  where  uic  =  (Ko/M')11*,  which  she 
prefers  to  the  second  order  theory.  With  the  zero  order  theory,  the  curves  shown  in  figures  3  and  14 
become  straight  lines  coincident  with  the  lines  in  the  figures  for  small  radii-  Trehu  points  out  that 
the  flexibility  of  the  OBS  base  can  lead  to  lower  damping  and  a  marked  decrease  in  the  frequency 
dependance  of  K.  Also,  anelastic  attenuation  within  the  sediments  should  lead  to  a  decrease  in  the 
variation  of  the  spring  and  damping  coefficients  with  frequency.  However,  the  second  ordeq  theory 
seems  to  fit  our  (more  numerous)  data  better  than  the  zero  order  theory. 

3) .  If  the  density  of  the  OBS  is  the  same  as  that  of  the  surrounding  medium,  then  the  OBS  is 
neutrally  bouyant  and  the  OBS  should  move  with  its  surroundings  for  seismic  inputs  regardless  of 
coupling  effects.  An  instrument  with  a  low  vertical  profile  and  density  equal  to  that  of  the  sediment 
should  move  with  the  sediment.  This  is  not  true  if  the  OBS  dimensions  are  larger  than  a  fraction 
of  a  wavelength,  because  the  velocity  of  a  seismic  wave  through  the  OBS  will  not  be  the  same  as 
the  velocity  in  its  surroundings,  and  the  OBS  will  distort  the  motion.  This  distortion  is  especially 
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important  for  OBSs  on  sediment  since,  for  example,  if  v,  —  10  m/s,  the  wavelength  at  20  Hz  is 
only  0.5  m.  Thus,  even  with  neutral  density,  distortion  can  be  expected  at  higher  frequencies. 

The  adverse  effects  of  coupling  increase  as  the  OBS  density  increases  above  the  sediment  density 
as  shown  in  Figure  7.  If  care  has  been  taken  to  insure  that  all  coupling  resonances  are  well  above 
the  frequencies  of  interest,  then  density  will  not  be  a  problem,  but  designing  for  neutral  density 
will  add  insurance  for  having  a  predictable  and  smooth  transfer  function. 

4) .  Low  profile  and  small  cross  section  in  the  water.  In  many  situations,  the  ocean  water  and 
the  sediments  have  a  discontinuity  of  horizontal  particle  motion  at  the  ocean  floor.  This  is  obvious 
in  the  case  of  significant  ocean  currents  [Duennebier  et.al.  (1981),Trehu  and  Solomon  (1981)],  but 
this  discontinuity  can  also  be  a  problem  for  seismic  inputs,  especially  shear  waves.  If  most  of  the 
volume  of  the  OBS  is  in  the  water,  the  instrument  will  tend  to  move  with  the  water  rather  than 
with  the  sediment.  The  water  and  sediment  move  together  for  vertically  incident  compressional 
energy,  but  for  other  angles  of  incidence  and  wave  types  they  do  not.  The  bottom  of  an  OBS  with  a 
large  cross  section  in  the  water  will  tend  to  move  with  the  sediment  while  the  top  will  tend  to  move 
with  the  water.  This  torquing  will  most  likely  result  in  detected  motion  that  is  severely  distorted, 
with  horizontal  particle  motion  recorded  by  the  vertical  sensors.  If  the  vertical  profile  of  the  OBS  is 
as  low  as  possible  and  the  base-to-height  ratio  is  large,  this  problem  should  be  eliminated.  A  lower 
vertical  profile  will  produce  a  lower  rocking  mode  sensitivity  and  a  higher  rocking  mode  coupling 
frequency.  However,  increasing  the  vertical  cross  section  within  the  sediment  should  improve  the 
response  to  horizontal  particle  motion. 

5) .  OBS  symmetry  and  minimum  density  gradients  will  tend  to  reduce  distortion  caused  by 
cross  coupling.  An  OBS  with  its  sensors  located  at  the  vertical  axis  of  symmetry  between  the 
center  of  buoyancy  and  the  center  of  mass  will  give  a  minimum  of  cross  coupling  problems.  As 
noted  earlier,  the  center  of  buoyancy  and  center  of  mass  should  be  as  dose  as  possible  to  avoid 
torquing  due  to  density  gradients  in  the  package.  If  the  OBS  is  elongated  in  one  horizontal  direction 
it  will  be  less  sensitive  to  cross  coupling  in  that  direction.  If  the  geophones  are  located  off  the  axis 
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of  symmetry  they  will  be  more  sensitive  to  rocking  and  twisting  modes  of  motion  of  the  package 
(Luco  (1974)).  Changes  in  shape  in  any  direction  will  modify  the  response,  as  virtual  mass  is  very 
shape  dependent,  thus  maximizing  symmetry  is  desirable. 

4.1  Practical  Considerations 

The  geophones  must  be  separated  from  the  necessarily  massive  package  that  contains  power,  elec¬ 
tronics,  recording,  and  recovery  equipment  to  obtain  OBS  masses  that  are  as  small  as  possible. 
Although  this  separation  may  add  some  complexity,  and  possibly  lower  reliability  than  when  ev¬ 
erything  is  in  one  package,  we  believe  that  the  advantages  of  the  increase  in  signal  fidelity  due  to 
minimization  of  coupling  problems,  the  likely  decrease  in  noise  from  water  currents  and  recording 
equipment  [Bycroft  (1956),Duennebier  et.al.  (1981)],  and  the  flexibility  in  sensor  configurations, 
far  outweigh  the  disadvantages. 

Separation  of  the  geophones  and  the  recording  package  requires  an  electro-mechanical  connec¬ 
tion  between  the  sensors  and  the  recording  package.  A  preamplifier  is  desireable  for  each  sensor 
in  the  geophone  package  to  prevent  problems  caused  by  electrical  leakage  to  sea  water  in  the  ca¬ 
ble  and  connectors.  The  cable  is  also  a  potential  problem  if  it  should  get  caught  under  the  OBS 
ballast,  and  thus  prevent  the  instrument’s  return.  This  problem  can  be  reduced  by  replacing  the 
traditional  solid  anchor  with  particulate  ballast  that  the  cable  can  pass  through  when  the  ballast  is 
released.  The  recording  package  itself  can  be  a  source  of  seismic  noise  and  signal  distortion  through 
reradiation  of  current  and  seismic  inputs.  Careful  design  and  adequate  separation  from  the  sensor 
package  can  eliminate  this  problem. 

As  geophone  packages  get  smaller,  other  problems  are  likely  to  be  encountered  because  of 
the  decreased  mass.  One  of  these  problems  is  the  electro-mechanical  connection  to  the  recording 
package.  This  link  must  be  carefully  designed  to  prevent  mechanical  forces  on  the  cable  from  being 
transmitted  to  the  sensor  package. 
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5  CONCLUSIONS 

y 

In  this  paper  we  have  presented  theoretical  arguments  and  data  that  suggest  that  ocean  bottom 
seismometers  should  be  designed  to  be  as  small  as  possible  to  prevent  poor  signed  fidelity  caused 
by  coupling  problems.  In  addition,  there  is  an  optimum  surface  area  in  contact  with  the  sediment 
depending  on  the  OBS  mass,  shape,  and  sediment  density.  It  is  found  that  the  frequency  above 
which  seismic  signals  are  distorted  is  directly  proportional  to  the  sediment  shear  velocity.  We 
conclude  that  significant  improvement  in  the  quality  of  seismic  data  can  be  made  if  the  guidelines 
presented  in  this  paper  are  followed  in  the  design  of  future  ocean  bottom  seismometers. 
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7  figure  captions 


Figure  1.  Fi  and  F%  vs.  dimensionless  frequency  (o)  from  equations  5  for  a  circular  disk  in  contact 
with  the  sediment.  The  points  on  the  curves  are  from  an  analytical  solution  by  Bycroft  (1956). 
The  optimum  dimensionless  frequency  is  shown  by  the  dashed  line. 

Figure  2.  Effective  mass  ( M'/Mi  =  and  damping  coeficient  (h  =  vs.  dimensionless 

resonant  frequency  ac  for  a  circular  disk-shaped  OBS.  The  optimum  coupling  frequency  occurs  at 
ac  =  1.16,  shown  by  the  dashed  line. 

Figure  3.  The  ratio  of  coupling  frequency  to  sediment  shear  velocity  vs.  bearing  radius  for 
circular  disks  with  different  masses  (Mi).  When  the  OBS  is  immersed  in  water,  M;  is  replaced 
by  the  dynamic  mass,  AT*  =  A/,  +  M„.  The  curves  terminate  at  the  limit  of  reliable  second  order 
theory. 

Figure  4.  Increase  in  and  decrease  in  fc  max  with  increasing  mass  for  a  circular  disk  from 
Equations  12,  and  13. 

Figure  5.  The  ratio  of  motion  of  a  sphere  of  density  pi  to  that  of  surrounding  water,  vw,  as 
a  function  of  the  ratio  of  sphere  mass,  Mi,  to  the  mass  of  water  displaced  by  the  sphere,  Af„,  from 
equation  14.  When  the  sphere  density  is  zero,  its  motion  is  three  times  that  of  the  water. 

Figure  6.  The  ratio  of  motion  of  a  disk,  Vi,  to  that  of  surrounding  water,  vw,  as  a  function  of 
the  ratio  of  disk  radius  r,  to  the  disk  thickness,  t,  from  equation  15.  The  water  is  moving  in  a 
direction  perpendicular  to  the  plane  of  the  disk. 

Figure  7.  Theoretical  seismic  phase  and  amplitude  response  for  three  values  of  a  for  h  =  0.1 
from  equations  21  and  22.  These  curves  are  do  not  depend  on  OBS  shape. 

Figure  8.  Theoretical  seismic  phase  and  amplitude  response  for  three  values  of  h  for  a  =  0.5 
from  equations  21  and  22. 

Figure  9.  Schematic  drawing  of  shaker  table  coupling  experiment.  The  table  motion  is  mon¬ 
itored  with  geophone  R,  and  the  motion  of  a  sphere  suspended  by  a  spring  is  monitored  by  an 
identical  geophone  (T)  inside  the  sphere.  A,  vertical  tests;  B,  horizontal  tests. 
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Figure  10.  Seismic  phase  and  amplitude  response  of  a  sphere  suspended  in  air.  The  smooth 
curves  are  solutions  to  Equations  21,  and  22  for  the  spheres  using  the  appropriate  parameters  from 
Table  1.  Figures  10a  and  b  are  for  vertical  motion  of  a  sphere  less  dense  than  water,  and  10c  and 
d  are  for  horizontal  motion  of  a  sphere  more  dense  than  water. 

Figure  11.  Seismic  phase  and  amplitude  response  of  a  sphere  suspended  in  water.  The  smooth 
curves  are  solutions  to  Equations  21  and  22  for  each  sphere  using  the  appropriate  parameters  from 
Table  1.  Figures  11a  and  b  are  for  vertical  seismic  motion  of  a  sphere  less  dense  than  water,  and 
11c  and  d  are  for  horizontal  seismic  motion  of  a  sphere  more  dense  than  water. 

Figure  12.  Example  of  the  comparison  between  theoretical  response  of  an  OBS  to  a  vertical 
transient  test,  and  the  actual  response  of  a  geophone  (dashed  line  is  from  Equation  24).  The  noise 
spectrum  is  shaded.  The  time  domain  curve  is  also  shown.  Care  must  be  taken  in  these  comparisons 
to  correct  for  the  geophone  response. 

Figure  13.  Theoretical  ratio  of  the  coupling  frequency  to  the  sediment  shear  velocity  (/c/t>,), 
compared  with  observed  coupling  frequencies,  (/c),  from  vertical  transient  tests.  fc/v,  is  calculated 
using  equation  6  and  appropriate  parameters  from  Table  2.  Each  symbol  represents  the  average 
of  several  vertical  transient  tests  for  a  particular  OBS  at  either  Lopez  Island  (squares),  the  1981 
Keehi  Lagoon  test  (triangles),  or  the  1982  Keehi  test  (diamonds).  The  symbol  size  is  proportional 
to  the  log  of  the  appropriate  dynamic  mass  from  Table  II.  Each  point  can  be  used  to  estimate  the 
shear  velocity  of  the  sediment  below  the  OBS.  The  two  lines  are  iso- velocity  lines  at  15  and  30  m/s. 

Figure  14.  /c//c  vs.  the  ratio  of  the  radius  of  the  area  in  contact  with  the  sediment 

to  the  optimum  radius.  The  data  points  are  from  Table  2  (assuming  that  p,  =  1570 kg/m3  and 
v,  =  20m/ s),  and  the  theoretical  curve  obtained  from  equations  6,  12  and  13,  is  independent  of 
sediment  parameters.  fc/fc  max  —  1.73[j^y^y]1/2.  Dashed  line  is  the  limit  of  reliable  second 
order  theory. 

Figure  15.  Theoretical  Q  from  Equation  7  versus  measured  Q.  Note  that  high  values  of  Q  are 
observed  only  from  large  OBSs. 
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ABSTRACT 


Ocean  bottom  seismographs  (OBS)  have  been  in  use  for  almost  50 
years.  The  earliest  instruments  were  used  for  short  range  explosion 
refraction  experiments  but  were  soon  abandoned  when  techniques  involving 
near  surface  sources  and  receivers  were  introduced.  Development  of  OBS 
technology  revived  during  the  late  1950's  inspired  by  interest  in  the 
detection  and  identification  of  clandestine  nuclear  explosions  and  sup¬ 
ported  largely  by  the  VELA-UNIFORM  program  of  the  United  States  govern¬ 
ment.  Since  that  time,  OBS  research  has  expanded  to  include  over  25 
organizations  in  several  countries.  Generally  each  group  has  developed 
its  own  instrumentation  and  experimental  procedures.  Instruments  fall 
into  four  classes:  (1)  self-contained,  free-fall,  pop-up;  (2)  acousti¬ 
cally  telemetering  to  near  surface  for  retransmission  or  recording;  (3) 
electrically  telemetering  via  cable  to  the  surface;  and  (4)  permanently 
linked  to  shore  by  cable.  Most  fall  in  class  (1).  Sensor  configura¬ 
tions  range  from  a  single  vertical  component  (or  hydrophone)  to  three- 
components  plus  hydrophone;  seismometer  free  period  is  generally  one 
second  or  less;  frequency  response  is  generally  limited  to  100  Hz  or 
less.  Instrument  masses  range  from  60  to  600  kg,  recording  is  analog 
(AM  or  FM)  continuous  or  digital  (time  windowed  and/or  triggered); 
dynamic  range,  data  storage  capacity  and  deployment  time  all  cover  a 
wide  range  of  values.  Most  applications  are  for  deteiminat ion  of  velo¬ 
city  structure  and/or  local  seismicity.  There  has  teen  lesser  interest 
in  teleseismic,  engineering  or  strong-motion  applications.  Major  prob¬ 
lem  areas  or  areas  where  improvements  are  needed  include  bottom-coupling 
(signal  distortion)  and  noise;  data  storage;  system  reliability,  espe¬ 
cially  over  long  deployment  times;  ease  of  handling  during  launch  and 
retrieval;  and  ease  of  shipboard  checkout  before  and  after  deployment. 
OBS  have  recently  been  supplemented  by  ocean  subbottom  seismometers 
emplaced  beneath  the  sediment  in  a  bore-hole. 


INTRODUCTION 

OBS  operated  at  selected  oceanic  sites  afford  an  opportunity  to 
explore  seismic  wave  propagation  under  the  deep  ocean  basins  and  across 
continental  and  island  margins;  ambient  seismic  noise  in  the  ocean 
basins,  both  natural  and  cultural;  the  geneiation  of  storm  microseisms; 
and  the  distribution  and  mechanisms  of  local  submarine  earthquakes.  In 
addition,  OBS  are  required  for  uniform  areal  obse-vations  of  seismic 
waves  from  teleseismic  events.  As  mentioned  in  the  abstract,  by  far  the 
majority  of  OBS  have  been  of  the  pop-up  variety,  designed  and  used  for 
relatively  short  term  active  and/or  passive  experiments. 


821 


The  design  of  a  good,  reliable  OBS  is  a  challenging  task.  In  addi¬ 
tion  to  the  requirements  of  any  remotely  operating  seismograph,  the 
instrument  must  be  designed  to  be  launched,  retrieved,  and  refurbished 
from  a  rolling,  often  wet  (with  salt  water)  ship;  it  must  couple  well  to 
a  (very!)  soft  or  hard  bottom  and  poorly  to  flowing  near-bottom  water. 
There  is  still  much  room  for  improvement  over  even  the  best  of  the 
existing  OBS.  As  is  the  case  for  land  instruments,  probably  more  so, 
optimum  design  is  a  strong  function  of  the  observational  objectives. 
For  example,  seismic  noise  on  the  ocean  floor  is  low  above  about  2  Hz 
and  often  very  high  in  the  storm  microseism  band,  three  to  10  seconds. 
It  is  still  uncertain  how  much  the  data  are  influenced  by  poor  bottom 
coupling.  Partially  because  of  the  high  noise  levels,  in  addition  to 
instrumental  difficulties,  relatively  little  work  has  been  done  with 
long-period  instruments  and  most  of  the  current  interest  is  above  1  Hz. 
However,  interest  in  the  long-period  range  is  growing.  At  the  other  end 
of  the  spectrum  there  is  interest  in  extending  the  frequency  band  to  a 
few  hundred  (or  thousand)  Hz  for  high  resolution  sediment  studies  and 
acoustic  studies  of  Navy  interest. 

Most  existing  OBS  are  designed  for  the  detection  of  small  signals, 
well  within  the  linear  range  of  bottom  motion,  and  only  such  instruments 
will  be  discussed  in  this  paper.  Steinmetz  et  al.  (1979)  describe  an 
OBS  specifically  designed  for  recording  strong  motion. 

An  excellent  review  paper  by  Whitmarsh  and  Lilwall  (1983)  provides 
a  fairly  complete  general  description  of  OBS  developed  and  used  through 

1981.  It  includes  an  extensive  bibliography  of  OBS  related  papers,  many 
of  which  provide  specification  tables  and  schematic  drawings  of  opera¬ 
tional  OBS.  A  more  recent  review  paper  by  Prothero  (1984)  discusses  the 
state  of  the  art  of  OBS  technology  through  1983,  based  partly  on  a 
workshop  held  at  the  University  of  California,  Santa  Barbara  in  July 

1982.  Since  1970,  a  sequence  of  papers  on  OBS  development  has  been  pub¬ 
lished  in  Marine  Geophysical  Researches.  Most  of  these  papers  are 
listed  by  Whitmarsh  and  Lilwall  (1983).  Physical  characteristics  of  the 
12  operational  OBS  deployed  during  the  Lopez  Island  OBS  intercomparison 
experiment  (conducted  in  June  1978)  are  tabulated  in  Sutton  et  al. 
(1981b).  Schematic  drawings  and  additional  specifications  for  these 
instruments  are  contained  in  Hawaii  Institute  of  Geophysics  Technical 
Report,  HIG-80-4,  Hawaii  Institute  of  Geophysics,  Honolulu. 

In  this  paper  I  shall  not  attempt  another  comprehensive  review  of 
the  field  but,  rather,  shall  limit  the  discussion  to  a  brief  review  of 
the  history  and  the  current  status  with  which  I  am  most  familiar, 
emphasizing  research  and  instrument  development  subsequent  to  that 
reported  in  the  Whitmarsh  and  Lilwall  and  the  Prothero  papers. 


EARLY  OBS 

Ewing  and  Ewing  (1961)  describe  the  history  and  status,  at  that 
time,  of  OBS  development.  The  earliest  OBS  wire  designed  for  small 
scale,  active  refraction  experiments  with  both  the  self-contained  OBS 
and  the  explosives  deployed  either  by  cable  o-  ti  e-fall  (pop-up).  In 
the  latter  case,  the  sensor  was  mechanically  isolated  : rora  the  flotation 
required  for  retrieval.  Figure  1  is  a  block  diagram  of  the  acous.ic 
telemetering  OBS  described  by  Ewing  and  Ewing.  Data  from  a  single  vert¬ 
ical  component  2  Hz  geophone  were  transmitted  by  frequency  modulation  of 
a  12  k  Hz  carrier  and  received  with  a  standard  echo-sounder  transducer. 
Different  mechanical  configurations  involved  either  cable  or  free-fall 
deployment  with  no  provision  for  retrieval  of  the  bottom  package.  Both 
explosion  and  earthquake  signals  were  well  recorded  with  this  system. 
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Figure  1.  System  block  diagram  of  acoustic  telemetering  OBS  (Ewing  and 
Ewing ,  1961). 
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Figure  2  shows  the  location  of  the  Columbia  (Pt.  Arena)  OBS,  about 
100  miles  west  of  San  Francisco,  California  in  four  kilometers  water 
depth,  that  operated  for  over  six  years,  1966-1972  (Sutton  et  al .  1965 
and  1969).  The  bottom  package  (Figure  3)  was  actually  a  small  geophysi¬ 
cal  observatory  containing:  three-component  long  (15  second)  and  short 
(one  second)  period  seismometers;  low  (crystal)  and  high  (coil)  fre¬ 
quency  hydrophones;  Vibratron  tidal  pressure  sensor;  crystal  temperature 
sensor;  and  water  current  magnitude  and  direction  sensors.  Ultra  long- 
period  (feed-back)  outputs  from  the  long-period  seismometers  provided 
tidal  gravity  and  tilt  data.  The  instrument  package  was  connected  to 
shore  at  Pt.  Arena,  California  via  two-conductor  deep-sea  cable.  Power 
and  30  command  channels  to  the  OBS  and  18  data  channels  from  the  OBS  to 
the  shore  station  were  transmitted  through  this  cable.  Fxcept  for  the 
mechanical  current  sensor,  which  failed  after  about  two  months  (probably 
from  biological  growth  and/or  silting)  all  instruments  operated  until 
the  cable  was  damaged  by  wave  activity  near  shore. 


Data  from  the  Columbia  OBS  are  included  in  a  number  of  reports  on 
short-  and  long-range  explosion  experiments  and  studies  of  local  seismi¬ 
city  and  teleseisms,  tides  and  tidal  currents,  and  microseisms.  One 
interesting  observation  is  that  the  OBS  often  recorded  offshore  earth¬ 
quakes  not  observed  by  land  stations  at  comparable  epicentral  distance 
and  the  converse  for  on-shore  earthquakes  (Sutton  et  al . ,  1969). 


Figure  4  illustrates  the  short-period  signals  from  a  local  earth¬ 
quake  recorded  by  the  Columbia  OBS  at  about  250  km  epicentral  distance. 
It  is  of  interest  to  note  that  each  of  the  three  motion  components  and 
the  pressure  component  emphasize  different  aspects  of  the  signals  pro¬ 
viding,  e.g.  information  on  direction  of  propagation  and  arrival  type. 
The  hydrophone  record  has  the  simplest  character;  a  strong  surface 
reflection,  P  ,  because  of  the  amplification  of  pressure  of  a  down¬ 
going  signal  (combined  with  its  bottom  reflection)  near  the  bottom;  and 
lack  of  the  shear  arrival,  Sn.  Pn  is  strong  on  the  vertical,  small  on 
(longitudinal  horizontal),  and  missing  on  Hj  (transverse  horizontal). 
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Figure  2.  Location  of  Columbia  OBS.  ARC,  UKI ,  and  BRK  are  stations  of 
the  Berkeley  (University  of  California)  seismic  network.  OBS  recording 
station  was  at  Pt.  Arena.  Mendocino  Escarpment  trends  E-W  near  40°N. 
Depth  contours  in  meters  (Sutton  et  al .  ,  1965). 


Figure  3.  Columbia  OBS  before  connection  to  cable.  Three-component  SP 
and  LP  seismometers  are  contained  in  two  of  white  spheres.  Other  sen¬ 
sors  are  mounted  on  frame.  Plastic  cover  on  current  meter  was  removed 
before  lowering.  Sphere  diameter  approximately  0.6  m  (Sutton  et  al., 


1965). 
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Figure  4.  P  and  S  arrivals  from  a  local  earthquake  recorded  by  the 
Columbia  OBS  (Figures  2  and  3).  Traces,  from  top,  are:  short-period 
hydrophone  (coil);  short-period  vertical  seismometer  (SPZ);  and  two 
orthogonal  short-period  horizontal  seismometers  ( SFH ^  and  SPH^).  Iden¬ 
tified  arrivals  are:  mantle  refracted  P  and  S  (Pn  and  Sn);  surface 
reflected  Pn  ( P  ) ;  and  Pn  converted  to  S  at  the  base  of  the  sedimentary 
layer  (Sp)  (Sutton  et  al.,  1969). 

AM  three  seismometers  record  the  converted  wave  Sp  and  the  shear 
arrival  Sn,  which  begins  earlier  on  the  vertical  and  longitudinal  hor¬ 
izontal  components  because  of  partial  conversion  to  P  Delow  the  ocean 
bottom.  The  resonant  character  of  the  seismometer  signals  results 
mainly  from  OBS-bottom  coupling  resonance,  to  be  discussed  more  fully 
later,  but  may  be  partly  due  to  sediment  velocity  structure.  Spectra 
from  P,  S,  and  T  phase  arrivals  on  a  given  seismometer  component  exhibit 
nearly  identical  resonant  peaks.  However,  the  resonant  frequency  is 
different  for  each  component:  7,  5,  and  6  Hz  for  Z,  H  and  respec¬ 
tively.  This  observation  strongly  suggests  that  tne  instruments  are 
recording  different  resonant  coupling  modes  excited  by  the  signals  (Sut¬ 
ton  et  al.,  1969).  Sutton  et  al.  (1965)  calculated  possible  resonant 
coupling  frequencies  for  the  Columbia  OBS  to  be  about  70  and  60  Hz  for 
vertical  and  rocking  motion,  respectively,  well  above  the  frequency  band 
of  interest  at  that  time.  That  calculation  was  based  upon  an  estimated 
sediment  shear  velocity  of  0.15  km/sec.  More  recent  observations  give 
shear  velocities  at  the  top  of  marine  sediments  an  order  of  magnitude 
lower.  Since  the  coupling  frequency  is  proportional  to  the  shear  velo¬ 
city,  a  velocity  of  15  m/sec  would  produce  calculated  frequencies  of  7 
and  6  Hz,  as  observed  in  the  data.  It  is  interesting  to  note  that  most 
of  the  operational  pop-up  OBS  tested  during  the  Lopez  Island  experiment, 
to  be  discussed  later,  exhibited  coupling  frequencies  between  about  5 
and  10  Hz  and  the  sediment  shear  velocity  near  the  bottom  was  measured 
to  be  15-20  m/sec. 

An  array  of  hydrophones  near  Wake  Island  and  connected  by  cable  to 
the  island  was  installed  more  than  25  years  ago  to  locate  nose  cone 
impacts  from  intercontinental  ballistic  missile  tests.  Eleven  of  these 
phones  are  still  operational.  Six  are  on  the  ocean  bottom  at  5.5  km 
depth  and  make  up  the  vertices  and  center  of  a  pentagon  about  40  km 
across;  five  are  located  at  three  different  sites  near  the  axis  of  the 
SOFAR  channel  (about  1  km  depth).  The  combined  array  spans  more  than 
300  km.  Signals  from  eight  of  these  hydrophones,  five  bottom,  and  three 
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SOFAR,  have  been  recorJed  digitally  by  the  Hawaii  Institute  of  Geophy¬ 
sics  on  a  continual  basis  since  September  1982.  The  hydrophone  signals 
are  recorded  through  low-noise  low-frequency  preamplifiers  and  are  used 
for  seismological  studies.  The  array  is  unique,  both  in  its  deep  ocean 
location  and  in  its  great  distance  from  other  seismological  stations. 

McCreery  et  al.  (1983)  examine  the  spectral  charac teris t ics  of  P 
waves  from  shallow  focus  earthquakes  and  nuclear  explosions  and  samples 
of  ocean-bottom  background  noise  recorded  on  the  Wake  array.  Ambient 
noise  levels  on  the  ocean  bottom  near  Wake  are  comparable  to  levels  at 
the  quietest  continental  sites  for  frequencies  between  three  and  13  Hz. 
At  lower  frequencies  the  noise  increases  rapidly  toward  the  "storm" 
microseism  maximum  near  0.2  Hz  where  the  level  is  above  that  for  con¬ 
tinental  sites. 

SHALLOW  WATER  OBS 

In  a  series  of  papers  Brocher  et  al.  (Brocher,  1983;  Brocher  et 
al.,  1981  and  1983)  reported  on  the  results  of  a  nine  day  OBS  experiment 
conducted  on  the  continental  margin  off  Nova  Scotia  in  June  1975.  Three 
HIG  telemetering  OBS  were  placed  in  a  linear  array  along  the  dip  of  the 
bottom  with  approximate  spacing  of  18  km  at  depths  of  67,  140,  and  L 30 1 
m.  The  lack  of  resonance  in  the  recorded  signals  suggests  that  the 
low-profile  instrument  packages  (7:1  length- to-d iameter  ratio)  provided 
excellent  coupling  to  the  sediments  for  the  three-component  2  Hz  natural 
frequency  geophones  and  avoided  the  coupling  of  the  seismometers  to 
ocean-bottom  currents.  The  OBS  were  hard-wired  directly  to  the  record¬ 
ing  ship  and  were  anchored  to  the  seafloor  by  a  heavy  communication 
cable.  The  broadband  frequency  response  of  both  the  incorporated  hydro¬ 
phone  and  the  geophones  together  with  the  FM  recording  of  the  signals  on 
(1.5  IPS)  tape  drives,  contributed  to  the  high  quality  of  the  recorded 
data.  The  horizontal  geophones  in  the  OBS  were  oriented  along  (horizon¬ 
tal  longitudinal)  and  perpendicular  (horizontal  transverse)  to  the  long¬ 
est  axis  of  the  OBS  package  (and  the  cable;. 

Data  from  a  variety  of  seismic  sources  were  recorded  during  the 
nine  day  invest igat ion .  The  purpose  of  the  experiment  was  to  study  the 
propagation  of  noise  produced  by  a  drilling  rig  platform  located  approx¬ 
imately  18  km  shoreward  from  the  shallowest  OBS.  Signals  and  background 
from  air  gun  profiling,  SUS  charges,  shipping,  biological  activity,  bot¬ 
tom  current,  varying  weather  conditions,  and  seismicity  were  investi¬ 
gated,  in  addition  to  the  rather  low  level  of  signal  observed  from  the 
drilling  operations.  P  and  S  velocity  structure  in  the  subbottom  was 
determined;  S/N  ratios  for  the  hydrophones  and  geophones  in  the  band  one 
t«>  30  Hz  were  compared;  and  low  frequency  propagation  along  strike  and 
up  and  down  dip  of  the  topography  was  studied.  Shear  velocities  were 
obtained  from  analyses  of  both  Stoneley  (Scholte)  and  refracted  shear 
waves . 


Related  studies  were  reported,  at  a  previous  NATO  conference  on 
bot tom- i nte rac t ing  ocean  acoustics,  by  Rauch  (1980)  and  Schirmer  (1980). 
In  both  of  these  studies  tethered  OBS  were  utilized  to  detect  interface 
waves  in  shallow  water. 

BOREHOLE  OBS 

A  number  of  experiments  have  recently  been  conducted  using  seismom¬ 
eters  clamped  in  a  drillhole  below  the  ocean  floor.  In  these  experi¬ 
ments  thf»  drilling  ship  Glomar  Challanger  of  the  Deep  Sea  Drilling  Pro¬ 
ject  was  used.  Two  papers  based  on  data  from  one  such  installation, 
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located  about  400  miles  southeast  of  the  central  Kuril  Islands  are 
included  in  this  volume  (Duennebier  et  al.,  1985;  and  Carter  et  al., 
1985).  This  ocean  subbottora  seismograph  (OSS)  is  emplaced  21  m  below 
the  sediment-basalt  boundary  at  a  depth  of  5.85  km.  The  OSS  includes  a 
self-contained  power  and  recording  system  that  can  be  left  on  the  ocean 
floor  and  retrieved  and  refurbished  at  a  later  time  (Duennebier  and 
Blackinton,  1983).  Stephen  (1981,  and  earlier  papers)  has  used  three- 
component  borehole  data  from  shots  recorded  on  the  drilling  ship  via 
direct  cable  connection  to  investigate  anisotropy  in  the  oceanic  crust. 
Carter  et  al.  (1984)  compared  data  from  an  OSS  emplaced  in  soft  sediment 
at  a  depth  below  bottom  of  194  m  with  that  from  an  HIG  ISOBS  (described 
later  in  the  section  on  Isolated  Sensor  OBS).  The  OSS  and  OBS  exhibited 
remarkably  similar  S/N  characteristics.  Absolute  noise  levels  were 
similar  above  10  Hz  but  below  10  Hz  the  OSS  was  quieter  by  three  to 
eight  dB. 

Another  borehole  seismograph,  the  Marine  Seismic  System  (MSS),  is 
similar  in  design  concept  but  more  elegant  that  the  OSS  (Ballard  et  al., 
1984).  Shearer  and  Orcutt  (1985)  de  tected  and  analyzed  anisotropy 
within  the  oceanic  lithosphere  using  MSS  and  OBS  data  collected  during 
the  NGENDEI  seismic  experiment  in  early  1983,  about  500  miles  east  of 
the  Kermadec  Trench  in  the  southwest  Pacific. 

Absolute  noise  levels  for  the  OSS  and  MSS,  when  clamped  in  com¬ 
petent  rock  beneath  the  sediment,  are  significantly  lower  than  for  OBS 
on  adjacent  ocean-bottom,  especially  for  the  horizontal  components  of 
motion.  S/N  also  appears  to  be  higher  for  most  transient  signals. 


THE  ROSE  EXPERIMENT 

The  Rivera  Ocean  Seismic  Experiment  (ROSE)  was  a  large  combined  sea 
and  land  seismic  program  utilizing  both  explosive  sources  and  earth¬ 
quakes  to  study  the  structure  and  evolution  of  a  mid-ocean  ridge,  a 
major  oceanic  fracture  zone  and  the  transition  region  between  ocean  and 
continent  (Ewing  and  Meyer,  1982).  Because  permission  was  not  received 
to  conduct  the  experiment  in  Mexican  territorial  waters  as  originally 
planned,  the  marine  portion  was  relocated  from  the  Rivera  Fracture  Zone 
to  the  East  Pacific  Rise  south  of  the  Orozco  Fracture  Zone.  The  two 
month  field  program  ( January-March  1979)  involved  many  cooperating 
research  groups  and  a  large  number  of  OBS.  Preliminary  results  are 
given  in  a  series  of  papers  in  the  same  number  of  Journal  of  Geophysical 
Research  as  the  Ewing  and  Meyer  overview  paper  (Project  ROSE  Special 
Issue,  J.  Geophys.  Res.,  87,  B10,  1982).  These  and  subsequent  papers 
(e.g.  Bratt  and  Purdy,  1984;  Bratt  and  Solomon,  1984;  Trehu  and  Purdy, 
1984)  provide  a  great  deal  of  information  on  the  P  and  S  crustal  velo¬ 
city  structure  of  the  Orozco  Fracture  Zone  and  East  Pacific  Rise  in  that 
vicinity  and  on  the  local  seismicity.  In  addition  to  the  scientific 
results,  ROSE  has  provided  useful  lessons  for  possible  future  large 
cooperative  efforts  and  mod  if icat ions  in  instrument  design  concerning 
such  disparate  items  as  data  archiving  and  distribution;  optimum  joint 
location  of  instruments  and  ships;  instrument  turn-around  time  and  reli¬ 
ability;  timed  versus  acoustic  recall;  and  OBS  data  capacity  and  timed 
versus  triggered  versus  continuous  recording. 

A  paper  by  Rowlett  and  Forsyth  (1984),  which  discusses  microearth¬ 
quakes  observed  at  the  intersection  of  the  Verna  Fracture  Zone  and  the 
Mid-Atlantic  Ridge,  provides  an  interesting  comparison  with  similar  stu¬ 
dies  from  ROSE  on  the  more  rapidly  spreading  East  Pacific  Rise. 
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ISOLATED  SENSOR  OBS 


In  most  of  the  OBS  used  today,  the  sensors,  electronics,  recorders, 
flotation,  and  ballast  are  contained  in  one  rigid  package.  This  confi¬ 
guration  requires  a  relatively  large  mass  and  large  vertical  cross  sec¬ 
tion  in  the  water  that  can  result  in  poor  coupling  to  the  bottom,  and 
increased  noise  level  from  response  to  water  flow  (Duennebier  et  al., 
1981).  The  use  of  an  isolated  sensor  can  minimize  these  problems.  In 
addition,  separation  between  the  sensors  and  the  main  OBS  package 
reduces  the  effect  of  possible  vibrations  from  the  data  recorder.  An 
isolated  sensor  OBS  (ISOBS)  developed  at  Hawaii  Institute  of  Geophysics 
(Byrne  et  al.,  1983)  is  shown  in  Figure  5.  The  cylinder  containing  the 
geophones  is  held  rigidly  to  the  main  package  during  free-fall  and  is 
released  after  arrival  at  the  bottom  by  the  solution  of  i  magnesium 
retaining  pin.  Direct  comparison  between  signals  recorded  by  a  single 
package  OBS  of  earlier  HIG  design  and  the  ISOBS,  deployed  near  each 
other  in  deep  water,  clearly  demonstrated  the  advantages  of  the  isolated 
sensor  configuration.  Additionally,  at  the  Lopez  Island  tests,  the  two 
OBS  with  separated  sensors  (ISOBS  and  another  supplied  by  MIT)  exhibited 
bottom  coupling  characteristics  superior  to  those  of  the  single  package 
OBS. 

Although  data  suggest  that  the  ISOBS  is  a  considerable  step  toward 
improving  the  fidelity  of  OBS  recording,  an  improved  prototype  isolated 
sensor  has  since  been  developed  at  HIC  (Figure  6)  known  as  the  "Opihi" 
(after  an  Hawaiian  limpet  that  adheres  to  underwater  surfaces).  The 
"Opihi  package  attempts  to  optimize  the  three  principal  design  parame¬ 
ters  found  to  affect  sensor  coupling  to  the  ocean  bottom  (Sutton  et  al., 
1981a, b).  The  parameters  are:  (1)  a  plate  with  a  relatively  large  sur¬ 
face  area  for  low  bearing  pressure,  (2)  nearly  neutral  density,  allowing 
the  package  to  "float"  in  the  sediments  and  (3)  a  buried  package,  for 
greater  coupling  to  solid  material  and  less  exposure  to  ocean  bottom 
currents , 

The  "Opihi"  package  deals  with  these  three  parameters  by  mounting 
the  ISOBS  geophone  pressure  case  on  a  plate  framework.  This  plate  is 
then  inserted  into  a  soft  mesh  bag  filled  with  a  loose  matrix  of  polyvi¬ 
nylchloride  pellets  (PVC)  or  pea  gravel.  This  method  distributes  the 
weight  of  the  pressure  case  over  a  wide  area,  making  it  about  neutral 
density  in  soft  sediments  (total  package  specific  gravity  of  about  1.4). 
The  loose  matrix  also  simulates  a  buried  sensor  and  should  appear  as  a 
small  bump  on  soft  sediments.  The  package  has  the  additional  advantage 
that  it  can  conform  to  bottom  contours  and  couple  well  on  hard  or  soft 
bottoms.  The  low  profile  of  the  package  makes  it  relatively  insensitive 
to  ocean  bottom  currents. 

In  essence,  this  sensor  package  should  have  good  coupling  and  low 
noise  generation  on  any  type  of  surface.  Initial  testing  of  the  "Opihi" 
sensor  package  design  has  shown  that  it  improves  coupling  over  the  pre¬ 
vious  ISOBS  sensor  package.  At  this  writing,  the  "Opihi"  has  not  been 
implemented  in  an  operational  OBS. 

CURRENT  AND  DESIRED  OBS  OPERATIONAL  SPECIFICATIONS 

In  1984,  Rondout  Associates,  Incorporated  conducted  a  survey  of  OBS 
developers  and  users.  Following  is  an  abbreviated  summary  Jeriveu  from 
the  34  responses  received  indicating  required  and  des  i  red  specifications 
(with  this  author's  comments  indicated  by  CHS): 

Major  Interests :  most  investigators  are  interested  in  Velocity 

Structure  and  Local  Seismicity.  Fewer  are  interested  in  Focal 
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Figure  3.  Isolated  sensor  OB S  geophone  release  scheme.  Separation  is 
approximately  one  meter  (Byrne  et  al.,  1983). 


Figure  6.  "Opihi"  isolated  sensor  package.  Scale  length  is  0.5  meters 
(  Byrne  et  al . ,  1983) . 


Mechanism  and  Exploration/Prospecting.  There  is  less  general 

interest  in  Teleseismic  Studies  and  Engineering  applications. 
Maximum  Depth :  most  want  capability  to  Normal  Ocean  (6000m)  depth; 
a  few  need  capability  to  work  in  Trenches  (12,000m);  fewer  still 
require  or  desire  Intermediate  (2000m)  or  Continental  Shelf  (500m) 
maximum  depth.  The  interest  in  shallower  depth  coincides  with 
Engineering  and  Strong  Motion  interests. 

Frequency  Range :  the  median  required  frequency  range  is  two  to  30 
Hz;  a  maximum  overall  range  of  0.01  to  100  Hz.  The  median  desired 
range  is  one  to  100  Hz  and  the  maximum  overall  range  is  DC-4  kHz ! 
Wider  band  instruments  emphasize  either  LP  or  HF,  typically  three 
decades . 

Components  Required :  most  require  or  desire  three-component 
seismometers  and  a  hydrophone.  Anything  less  is  generally  a 
compromise  for  data  storage  or  other  reasons.  (Each  of  the  four 
components  provides  unique  information.  However,  since  the  hydro¬ 
phone  generally  produces  relatively  clear  and  simple  records  with 
good  S/N,  relatively  undisturbed  by  poor  bottom  coupling  and 
current  noise,  it  is  sometimes  preferred  for  exploration  and  acous¬ 
tic  work.-GHS) 

Record ing :  everybody  requires  or  desires  digital  recording.  Some 
stick  to  analog  recording  (for  good  reason!-GHS)  because  of  the 
great  advantages  of  continuous  recording  in  the  marine  environment; 
limitations  of  digital  data  quantity;  problems  with  triggering  on 
noise,  and  not-tr igger ing  on  weak  signals;  and  logistical  problems 
with  timed  recording. 

Recording  Mode :  most  investigators  require  or  desire  a 
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Microprocessor  Controlled  system  and  prefer  the  Triggered  Mode  to 
the  Timed  Mode.  (Peal  and  Kirk,  1983,  describe  a  versatile 
microprocessor-controlled  OBS  data  logging  system.  This  subject  is 
also  covered  by  Prothero,  1984. -GHS)  A  significant  number  require 
(or  currently  have)  Continuous  Recording. 

Recording  Medium:  the  majority  currently  use  Cassette  Recorders; 
some  use  or  would  prefer  Reel-to-Reel  or  Cartridge  recorders.  Most 
of  those  using  digital  recording  would  desire  greater  capacity  than 
their  existing  system.  (The  recorder  is  a  major  problem  area.  In 
addition  to  problems  with  storage  capacity,  it  is  the  major  power 
drain,  is  one  of  the  least  reliable  components,  and  can  be  an 
important  source  of  vibration  noise.  The  OBS  community  would  wel¬ 
come  a  real  breakthrough  in  this  area. -GHS) 

Dynamic  Range:  the  median  required  dynamic  range  is  72  dB;  the 
responses  range  from  26  dB  (analog)  to  96  dB  without  gain  ranging 
and  130  dB  with  gain  ranging.  The  median  desired  dynamic  range  is 
120  dB;  the  responses  range  from  60  to  80  dB  and  40  dB  plus  AGC  to 
130  dB,  16  bit  floating  point. 

Deployment  Time:  the  median  required  Deployment  Time  is  30  days 
with  a  total  range  from  one  day  to  one  year.  The  median  desired 
deployment  time  is  six  months  with  a  range  of  one  day  to  "semi¬ 
permanent".  The  short  times  are  from  groups  principally  interesLed 
in  Shelf  and  Intermediate  Depth  Engineering  studies  and/or 
Exploration/Prospecting;  the  long  times  coincide  with  interest  in 
Teleseismic  and  Long  Period  work. 

Recording  T ime :  the  answers  to  this  item  are  given  in  terms  of  con¬ 
tinuous  recording  times  (analog),  the  number  of  short  event  time 
windows;  and/or  the  total  number  of  megabytes  (mB)  data  storage 
capacity.  The  median  values  listed  are  approximate.  The  median 
required  Recording  Capacity  is  about  20  mB,  digital;  or  20  to  30 
days,  analog.  The  required  values  range  from  20  minutes  (5  mB)  to 
90  mB,  digital;  and  four  Jays  to  two  months,  analog.  The  median 
desiied  Recording  Capacity  is  about  100  mB  digital  and  three  to  six 
months,  analog.  The  values  range  from  about  40  mB  to  oo  ,  digital, 
and  one  month  to  one  year,  analog.  (As  noted  before,  under  Record¬ 
ing  Medium,  this  is  an  area  that  needs  improvement . -GHS ) 

Acceptable  Clock  Drift :  (the  important  thing  is  uncorrectable 
drift.  Linear  drift  can  be  corrected. -GHS)  The  acceptable  uncer¬ 
tainty  depends  upon  the  experimental  objectives.  Total  uncertainty 
of  +  /-  0.1  second  or  less  is  required  for  most  earthquake  work. 

Short-period  array  studies  and  high  resolution  sediment  studies 
require  near  1  msec  accuracy.  Total  clock  stability  of  one  part  in 
10*  is  desirable. 

Maximum  Weight :  the  median  required  Maximum  Weight  is  150  kg  with 
values  ranging  from  45  kg.  to<  450  kg.  The  median  desired  maximum 
weight  is < 90  kg.  with  values  ranging  from  10  kg  to  300  kg..  A  few 
investigators  who  work  from  large  ships  are  relatively  unconcerned 
about  weight.  Workers  with  more  limited  ship  capabilities  and  per¬ 
sonnel  and/or  who  utilize  a  larger  number  of  instruments  are  con¬ 
cerned  to  keep  weight  and  handling  difficulties  to  a  minimum. 
Release :  everyone  prefers  an  acoustic  release,  generally  with  a 
timed  release  backup.  (Use  of  acoustic  telemetry  for  diagnostics 
is  useful.  Au  associated  transponder  used  for  precision  navigation 
is  invaluable  for  certain  experiments .-GHS) 

Bottom  Coupling:  (I  have  viewed  this  as  a  major  problem  area  for  a 
number  of  years. -GHS)  All  require  good  coupling  to  soft  bottoms  and 
most  to  rocky  bottoms.  Fewer  are  concerned  with  High  Currents. 
(Refer  to  further  discussion  below  under  Packaging 
Configurations .-GHS) 

Azimuthal  Orientation:  most  require  or  desire  an  independent  meas¬ 
urement  of  azimuth.  The  median  preferred  accuracy  is  +/-  5°  with  a 
range  of  values  from  0.5  to  10°. 
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Packaging  Configuration:  a  majority  of  the  respondents  require  or 
desire  either  Separated  Sensors  or  Low  Profile  or  both;  two,  indi¬ 
cate  (strong)  preference  for  a  single  package;  one  suggests  an 
adaptable  system  depending  upon  bottom  conditions.  (The  Lopez 
Island  experiment  indicated  that  the  coupling  performance  of  iso¬ 
lated  sensor  package  OBS  is  better  than  single  package  instruments, 
although  some  of  the  latter  performed  almost  as  well.  The  isolated 
sensor  instruments  produced  records  comparable  to  the  three  special 
"standard"  instruments  useJ  at  Lopez  Island  (Sutton  et  al.,  1981b). 

Use  of  a  well-designed  separated  sensor  package  generally  should 
improve  the  bottom  coupling  and  isolation  from  current  induced  and 
tape  recorder  noise.  It  should  also  be  a  more  flexible  design  per¬ 
mitting  various  sensor  packages  and  a  common  support  and  recording 
package.  The  requirement  of  including  flotation  and  reLLieval  aids 
on  a  single  package  instrument  makes  optimum  coupling  design  diffi¬ 
cult.  However,  generally  an  isolated  sensor  design  is  more  complex 
and,  therefore,  less  rel iable ,-GHS) 

Retrieval :  everyone  wants  a  Flasher  and  Radio  (although  the  Radio 
antenna  can  be  a  serious  source  of  current-induced  noise-GH3) . 
Several  require  or  desire  an  acoustic  pinger  and/or  transponder  to 
aid  in  retrieval. 

Maximum  Unit  Cost:  the  majority  listed  $10-20,000;  a  significant 
number  prefer  <$10,000;  or  will  accept  $20-40,000;  and  a  few  will 
go  to  >$40,000.  Generally,  the  higher  values  coincide  with  more 
stringent  OBS  requirements. 


Cenera 1  Evaluation  of  Current  OBS  Designs-(GHS) 

Major  problem  areas  or  areas  where  improvements  are  needed  include 
bottom-coupling  (signal  distortion)  and  noise;  data  storage;  system 
reliability,  especially  over  long  deployment  times;  ease  of  handling 
during  launch  and  retrieval;  and  ease  of  shipbound  checkout  before  and 
after  deployment.  Each  of  the  OBS  with  which  I  have  personal  familiar¬ 
ity,  principally  those  used  in  the  Lopez  and  ROSE  experiments,  has  some 
strong  and  weak  points.  The  external  configurations  of  many  of  the  OBS 
were  modified  subsequent  to  Lopez  to  improve  coupling  characteristics. 

The  trade  offs  between  analog  and  digital  have  been  discussed  ear¬ 
lier.  Most  investigators  prefer  the  digital  route.  The  limited  dynamic 
range  of  analog  recording  can  be  alleviated  by  gain  ranging  on  back¬ 
ground  noise,  multiple  gain  channels,  and  temporary  digital  storage  and 
subsequent  gain  range  before  analog  recording.  Much  digital  data  have 
been  lost  by  poor  trigger  algorithms  and  other  problems. 

Reliability  has  been  a  serious  problem  with  some  systems,  both  from 
lost  data  and  from  lost  instruments.  The  acceptable  loss  rate  depends 
upon  the  cost  and  uniqueness  of  the  data  and  the  instrument  package. 
Some  groups  have  experienced  very  low  rates;  and  others  10%  or  more. 
Occas  onally,  losses  come  in  batches.  Details  of  each  component  and  the 
system  as  a  whole  must  be  considered  carefully.  External  sensors 
(seismometers  and  hydrophones)  require  piessure  tight  electrical  lead- 
throughs;  cables  to  isolated  sensors  can  become  caught  cn  the  anchor  or 
bottom  rocks.  Release  mechanisms  can  fail;  glass  instrument  and  flota¬ 
tion  spheres  can  fail  after  repeated  use. 

Some  systems  are  difficult  to  launch  and  retrieve  without  damage; 
some  of  the  larger  instruments  require  special  handling  gear.  Some  sys¬ 
tems  require  extensive,  time-consuming  pre-launch  checkout,  limiting  the 
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number  of  instruments  that  can  be  deployed  for  a  given  experiment.  The 
conflicting  technical  requirements  or  desires  of  the  OBS  community  for 
low-power  and  low- ins t rument  volume,  versus  high  data  storage,  high 
dynamic  range  and  high  clock  stability  are  pushing  the  state  of  the  art. 


SIGNAL  DISTORTION  AND  NOISE  IN  OBS 

A  badly  designed  OBS  can  be  subject  to  severe  signal  distortion 
because  of  poor  coupling  to  the  soft  ocean  bottom  and  suffer  from  nigh 
levels  of  noise  caused  by  near-bottom  water  currents  shaking  the  instru¬ 
ment  frame.  Most  existing  operational  OBS  exhibit  some  amount  of  these 
undesirable  characteristics.  The  seismic  frequency  response  of  an  OBS 
frame  of  a  given  mechanical  configuration  coupled  to  sediment  of  a  given 
shear  velocity  (V^)  and  density  can  be  predicted  using  relationships 
that  have  been  tested  experimentally  using  in  situ  mechanical  tran¬ 
sients.  These  same  relationships  can  be  used  to  establish  guidelines 
for  designing  instruments  with  improved  coupling  response.  Seismic  sig¬ 
nals  couple  to  an  OBS  both  through  the  area(s)  of  contact  with  the  sedi¬ 
ment  and  directly  through  the  contact  with  the  water.  Assuming  that  the 
desired  response  is  to  faithfully  follow  the  motion  of  the  sediment  sur¬ 
face,  direct  coupling  to  the  water  for  hor izont a  1- trans l at  ion  and  rock¬ 
ing  modes  of  motion  should  be  minimized.  The  OBS  will  follow  the  sedi¬ 
ment  motion  for  frequencies  lower  than  the  natural  frequencies  of  the 
modes  of  coupling  motion.  Thus,  it  is  generally  desired  to  have  these 
coupling  frequencies  ( f  )  as  high  as  possible.  For  a  given  OBS  confi¬ 
guration,  there  is  an  optimum  effective  radius  for  the  area  of  coupling 
to  the  sediment  that  will  produce  the  highest  predictable  coupling  fre¬ 
quency.  The  optimum  radius  is  proportional  to  the  cube  root  of  the 
instrument  mass.  Using  these  optimum  radii,  the  coupling  frequency 
decreases  as  the  cube  root  of  the  mass,  i.e.  smaller  instruments  will 
have  higher  (more  favorable)  coupling  frequencies. 

Recognizing  the  importance  of  the  OBS-bottom  coupling  problem,  the 
Lopez  Island  OBS  experiment  was  conducted  in  Shoal  Bay,  Lopez  Island  in 
June  1978,  as  a  preliminary  to  ROSE,  in  order  to  compare  the  responses 
of  the  operational  OBS,  which  have  quite  varied  mechanical  configura¬ 
tions,  when  coupled  to  the  earth  through  soft  sediments  comparable  to 
those  of  the  ocean  floor,  and  to  determine  their  susceptibility  to  noise 
induced  by  near-bottom  currents  (Sutton  et  al.,  1981b).  (More  recently, 
a  similar  experiment  involving  a  different  set  of  OBS  was  conducted  in 
France;  I  am  not  aware  of  published  results.)  The  maximum  current  meas¬ 
ured  during  the  test  was  6  cm/sec  and  no  significant  correlation  was 
found  with  noise  on  the  OBS  during  the  experiment. 

Signals  recorded  on  OBS  geophones  are  often  narrow  band  and  pro¬ 
longed  (Figure  4).  Whether  these  features  are  natural  in  origin  or 
caused  by  the  presence  of  the  OBS  was  a  major  question  to  be  answered  by 
the  test.  Some  data  indicated  that  the  presence  of  the  OBS  on  the  soft 
sediments  could  severely  distort  the  motion  of  the  ground;  however,  the 
extent  of  this  problem  for  OBS  of  different  configurations  was  unknown 
The  complications  present  on  many  geophone  records  are  not  observed  on 
ocean  bottom  hydrophones,  thus  the  question;  do  geophone  data  adu  any 
useful  information  that  cannot  be  obtained  from  hydrophones? 

Some  data  had  indicated  that  near-bottom  currents  can  increase 
noi3e  levels  on  OBS  enough  to  make  them  unusable  during  high-current 
periods;  therefore,  we  wanted  to  test  as  many  systems  as  possible  for 
susceptibility  to  ocean  current  noise. 
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In  addition  to  12  operational  OBS  supplied  by  10  different  research 
organ i zat  ions ,  three  sets  of  three-component  (T  ®ls)  seismometers  and  a 
hydrophone  were  included  to  provide  standards  for  comparison:  a  "spike 
standard"  was  pushed  firmly  into  the  bottom;  a  "plate  standard"  with  a 
smooth,  hemispherical  superstructure  on  a  large,  flat,  circular  plate 
was  placed  on  the  sediments;  and  a  "neutral  density  standards"  with  a 
roughly  sphericaL  shape  was  floated  within  the  uppermost  sediment  The 
instruments  were  placed  within  a  few  meters  of  each  other,  and  four 
current  meters  were  installed  around  the  array  to  monitor  water  circula¬ 
tion  in  the  bay. 

Data  were  obtained  from  each  sensor  by  hard  wiring  analog  signals 
to  land  where  they  were  recorded  by  a  digital  system.  Seismic  records 
were  obtained  from  mechanical  transient  tests,  airgun  shots  and  blasting 
caps,  and  samples  of  background  noise.  Mechanical  transient  tests  were 
conducted  to  provide  some  quantitative  estimate  of  the  coupling  function 
between  the  ocean  bottom  and  the  OBS  package,  including  possible  cross 
coupling  between  horizontal  and  vertical  motions.  The  procedure  is 
analogous  to  the  classical  weight  lift  test  for  earthquake  seismographs, 
although  because  of  buoyancy  the  response  of  the  seismic  system  to  the 
mechanical  transient  test  is  not  the  same  as  that  to  a  seismic  input 
(button  et  al.,  1981a). 

An  Hawaii  Institute  of  Geophysics  Technical  Report  (HIG-80-4)  con¬ 
tains  a  large  number  of  records  from  each  recorded  sensor.  The  techni¬ 
cal  report  also  shows  the  designs  of  all  instruments.  All  the  data  of 
this  experiment  are  on  digital  tapes  and  may  be  requested  from  the  ROSE 
data  archive  at  Hawaii  Institute  of  Geophysics.  A  catalog  of  the  data 
appears  in  the  technical  report. 

Subsequent  to  the  Lopez  Island  experiment,  a  number  of  papers  on 
bottom  coupling  have  improved  our  understanding  of  the  problem  and  pro¬ 
vided  guidelines  for  better  OBS  design  (Sutton  et  al.,  1981  a  and  b; 
Zelikovitz  and  Prothero,  1981;  Lewis  and  Tuthill,  1981;  Trehu,  1985; 
Garmany,  1984;  Sutton  et  al.,  1985;  and  Sutton  and  Duennebier,  1985). 

Table  l  lists  parameters  that  must  be  considered  and  Table  II  sum¬ 
marizes  possible  problems  of  signal  distortion  and  noise  related  to 
OBS-bottom  coupling.  A  schematic  diagram  of  the  dynamic  coupling 
between  an  OBS  and  a  soft  bottom  for  vertical  motion  is  shown  in  Figure 
7.  In  general,  the  spring  "constant",  K,  is  not  equal  to  the  static 
value;  both  K  and  the  damping  "constant",  D,  are  functions  of  frequency. 
The  inertial  effect  of  the  water,  Mw* ,  and  sediment,  Ms*,  involved  in 
the  differential  motion  between  the  OBS  and  its  surroundings  may  also  be 
considered  to  be  frequency  dependent  (or  alternatively  included  in  the 
variation  in  K).  The  M*,  K,  and  D  are  also  shape  dependent. 

The  case  for  horizontal  input  is  oiiuilar  to  that  for  vertical 
motion,  with  the  added  complications  that:  there  can  be  two  independent 
inputs  in  each  orthogonal  direction,  one  through  the  water  and  one 
through  the  sediment;  and  cross-coupleu  rocking  motion  can  contribute  a 
major  portion  of  the  output  from  the  horizontal  seismometers  (and  from 
the  vertical,  if  it  is  located  far  from  the  center  of  rotation).  The 
response  of  a  horizontal  geophone  to  a  rocking  mode  tilt,  0,  is 
equivalent  to  a  horizontal  acceleration,  yt=gB,  so  that  the  total 
apparent  input  is  y+y^ . 

The  response  of  a  rigid  sphere  in  a  fluid  to  a  seismic  signal  in 
the  fluid  of  wavelength  large  compared  to  the  sphere  diameter  is  U/V- 
3/ ( 2 ( Mi /Mw) +1 )  where  U  and  V  are  the  motion  of  the  sphere  and  the 
incident  wave  in  the  fluid,  respectively,  and  Mi  and  Mw  are  the  mass  of 


TABLE  I 

OBS-BOTTOM  COUPLING  PARAMETERS 


A.  OBS 

1 .  Mass 

2.  Moments  of  Inertia 

3.  Structural  Rigidity 

4.  Symmetry 

5.  Shape/Smoothness 

6.  Cross-S(  't ion  to  Water 


7.  Cross-Section  to  Sediment 

8.  Density 

9.  Bearing  Area 

10.  Bearing  Pressure 
1 1  .  Burial 


SEDIMENT/ROCK 

1  . 

Density 

4. 

Heterogeneity 

2. 

3. 

Shear  Velocity 

Anelast icity 

5. 

Permeability 

WATER 

1  . 

2. 

Density 

V  iscos i ty 

3. 

Sound  Velocity 

D.  INTERFACES/STRUCTURE 

1.  Water-Sediment  Interface  Roughness 

2.  Sediment-Rock  Interface  Roughness 

3.  Velocity-Depth  Dependence 


TABLE  II 

OBS-BOTTOM  COUPLING  PROBLEMS 


SIGNAL  DISTORTION 


NOISE 


OBS  on  Sed  iment 

-Vertical  Output  from  Vertical 
Signal 

-Horizontal  Output  from  Horizontal 
Signal 

-Vertical  Output  from  Horizontal 
Signal 

-Horizontal  Output  from  Vertical 
Signal 

OBS  on  Rock 

-Rat t 1 ing 

-Large  Tilt 

-"Floating"  Footpad 

Direct  OBS-Water  Coupl ing 

-Discontinuous  Horizontal  Signal 
Across  Water/Sediment  Boundary 

-Inertial  Effect  (Most  OBS  have 
Large  Internal  Vertical 
Gradient  in  Density) 

Lateral  Heterogeneity  of  Bottom 

Tilts  on  Horizontals  from  Ultra-Low 
Phase  Velocity  Signals 


OBS  wn  Sediment 
-Amplification  of  Microseisms 
at  Bottom  Interface 
-Signal  Induced  "Noise"  from 
Irregular  Boundary 

-Water/Sediment 
-Sediment/Layer  2 


Direct  OBS-Water  Coupling 
-Current  Induced  Noise 
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Figure  7.  Schematic  of  OBS  coupling  to  ocean  bottom  for  vertical 
motion.  The  dashed  line  is  assumed  to  move  with  the  ocean  floor  adja¬ 
cent  to  the  OBS  (Sutton  et  al.,  l9Sia), 


the  sphere  and  the  displaced  water,  respectively.  Similar  relationships 
hold  for  bodies  of  other  shapes.  Since  current  OBS  designs  generally 
include  most  of  the  body  of  the  instrument  within  the  water,  direct 
seismic  coupling  between  water  and  OBS  is  certain  and  in  many  cases 
could  be  more  important  than  coupling  through  the  bottom  sediments.  From 
I  the  equation,  we  see  that  amplification  can  approach  a  factor  of  three 

for  very  low  sphere  densities  and  that  attenuation  occurs  for  densities 
greater  than  water.  As  expected,  when  Mi=Mw,  U=V.  Most  OBS  contain 
low-density  flotation  for  retrieval  in  their  upper  portion  and  more 
dense  disposable  anchoring  material  near  the  bottom-coupling  interface. 
This  density  gradienc,  in  addition  to  producing  torques  from  equal 
amplitude  horizontal  inputs,  will  modify  the  response  from  waves  having 
|  different  horizontal  amplitudes  in  the  water  and  in  the  sediment.  t'e.g. 

horizontally  propagating  compi e»» loud  1  waves  or  horizontally  polarized 
shear  waves).  The  vertical  component  is  less  complicated  since  signal 
amplitude  is  continuous  across  the  water-sediment  boundary  and  most  OBS 
have  appreciable  symmetry  about  a  vertical  axis.  However,  lateral  vari¬ 
ations  within  the  bottom  under  the  OBS  could  affect  this  component, 
also. 


A  simple  calculation  indicates  the  possibility  of  significant  noise 
resulting  from  response  to  tilts  produced  by  bottom  cuirents.  Consider 
a  spherical  OBS  with  cross  section  lm  ,  centered  in  the  water  1/2  m 
above  two  circular  footpads  of  radius  0.2  m  that  are  separated  by  1  m  in 
the  direction  of  the  currenjj,  resting  on  sediment  with  shear  velocity  20 
m/sec  and  density  1  g/cm  .  The  vertical  displacement  of  the  footpads 
(one  up,  one  down)^from  a  bottom  current  of  5  cm/sec  is  about  10  pm. 
This  is  roughly  10J  larger  than  the  vertical  signal  from  a  teleseismic  P 
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Figure  10.  OBS  resonant  coupling  frequency,  fc,  as  a  function  of  radius 
of  bot tom-coupl ing  area  and  instrument  mass.  Normalized  by  shear  velo¬ 
city  of  bottom,  Vg.  Curves  terminate  at  limit  of  second-order  theory 
(Sutton  and  Duennebier,  1985). 


wave.  The  horizontal  acceleration  associated  with  the  tilt  (y=g0)  is 
equivalent  to  a  horizontal  amplitude  of  200  m  at  6  second  period.  This 
value  is  about  20  times  larger  than  the  horizontal  signal  from  average 
"storm"  mic rose ismic .  Small  perturbations  of  bottom  currents  around 
such  large  mean  values  should  produce  significant  noise  on  most  opera- 

|  tional  OBS. 

Figures  8  and  9  show  the  amplitude  response  of  an  OBS  to  a  seismic 
input,  assuming  frequency  independent  coefficients,  for  two  values  of 
damping  and  three  values  of  the  coupling  coefficient  c.  It  can  be  seen 
that  the  OBS  witl  faithfully  follow  the  input  for  frequencies  well  below 
the  coupling  resonance,  u>c,  but  that  serious  distortion  can  occur  near 
resonance  and  significant  attenuation  can  occur  at  higher  frequencies. 

®  The  coupling  coefficient,  c,  depends  upon  the  OBS  density  (as  indicated 

earlier  for  a  sphere  in  water).  c=0  is  neutral  density;  and  gives  per¬ 
fect  coupling  independent  of  uic.  These  response  curves  are  reliable  to 
frequencies  somewhat  above  oj  c ,  but  there  is  some  question  regarding 
their  precision  at  higher  frequencies  (Trehu,  1985;  Sutton  and  Duenne¬ 
bier,  1985). 

|  The  relationship  from  second-order  theory  between  coupling  fre¬ 

quency  and  bearing  radius  (assuming  a  circular  bearing  area)  for  OBS  of 
different  mass  is  shown  in  Figure  10.  For  a  given  mass,  the  coupling 
frequency  increases  with  radius  to  a  maximum  value  near  the  limit  of 
validity  of  the  theory.  The  damping  also  increases  with  radius  and  is 
near  critical  at  the  maximum  of  these  curves.  Thus,  it  appears  that  the 
radius  of  tile  maximum  is  also  the  optimum  design  radius  for  a  given  OBS. 
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Note,  also,  that  the  larger  the  08S  the  more  difficult  it  is  to  obtain  a 
high  coupling  frequency. 


From  the  foregoing,  Sutton  and  Duennebier  (1985)  arrived  at  the 
following  desired  parameters  for  OBS  design:  (1)  minimum  mass,  (2) 
optimum  radius  in  contact  with  ocean  floor  for  a  given  instrument  mass, 
shape  and  sediment  density,  (3)  density  approximately  that  of  the  sedi¬ 
ment,  (4)  sensors  near  the  center  of  mass/vertical  symmetry,  and  (5) 
small  cross  section  within  the  water.  The  "Opihi"  (Figure  6)  satisfies 
these  requirements  quite  well. 

Garmany  (1984)  presents  a  generalized  procedure  for  recovering  true 
ocean  bottom  particle  motion  using  only  the  output  signals  from  a 
three-component  OBS  that  have  been  distorted  by  coupling  and/or  the 
local  bottom  conditions.  Sutton  et  al.  (1985)  point  out  some  problems 
with  Garmany's  analysis.  A  modified  version  of  Garmany's  method  might 
be  used  to  improve  the  quality  of  existing  OBS  data,  but  I  believe  it  is 
not  a  substitute  for  good  OBS  design. 
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Comment  on  "The  Recovery  of  True  Particle  Motion  From  Three-Component 
Ocean  Bottom  Seismometer  Data”  by  Jan  D.  Garmany 
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Chirniiiny  (1984J  presents  a  procedure  for  recovering  the 
true  ocean  bottom  particle  motion  using  the  output  signals 
from  a  three -component  ocean  bottom  seismometer  lOBSl.  In 
the  abstract  he  states.  "The  most  significant  implication  of  this 
work  is  that  accurate  ground  motion  can  be  recovered  from 
almost  any  three-component  OHS  ”  This  is  a  worthy  objective, 
since  it  is  well  known  that  much  OBS  data  are  degraded 
because  of  poor  coupling  through  soft  bottom  sediments  and 
water  The  paper  provides  an  interesting  generalized  approach 
that  might  be  used  to  improve  OBS  data  if  some  serious 
errors  and  omissions  contained  in  the  paper  arc  properly  cor¬ 
rected  and  considered 

tuirmuni  [1984J  assumes  that  the  OBS  has  three  inputs  and 
three  outputs  and  that  the  OBS  site  system  is  an  all-poJe 
minimum  delay  filter  He  indicates  that  these  two  assumptions 
are  required  in  order  lor  the  procedures  outlined  in  the  paper 
|  to  be  valid.  However,  in  general,  an  OBS  has  five  independent 

rectilinear  seismic  inputs:  one  vertical  input  and  four  horizon¬ 
tal  inputs,  two  from  the  bottom  sediment  and  two  from  direct 
coupling  to  the  near-bottom  water,  which  moves  differently 
from  the  bottom  for  various  seismic  wave  types.  Additionally, 
the  transfer  function  for  a  seismic  input  contains  one  or  two 
zeros  as  well  as  poles.  We  now  consider  further  each  of  Gar- 
9  many's  two  assumptions. 

The  response  of  a  rigid  sphere  in  a  fluid  to  a  seismic  signal 
in  the  fluid  of  a  wavelength  which  is  large  compared  to  the 
sphere  diameter  is  i  f  -  }  t2l.\f,  *■  I).  where  t  and  V 

are  the  motion  of  the  sphere  and  the  incident  wave  in  the 
fluid,  respectively,  and  M,  and  Af„  are  the  mass  of  the  sphere 
and  the  displaced  water,  respectively  [ Batchelor.  1967].  Simi- 
*  lar  relationships  hold  for  bodies  of  other  shapes.  Since  current 

OBS  designs  generally  include  most  of  the  body  of  the  instru¬ 
ment  within  the  water,  direct  seismic  coupling  between  water 
and  the  OBS  is  certain  and  in  many  cases  could  be  more 
important  than  coupling  through  the  bottom  sediments.  From 
the  equation  we  see  that  amplification  can  approach  a  factor 
|  of  3  for  very  low  sphere  densities  and  that  attenuation  occurs 

for  densities  greater  than  water.  As  expected,  when  Af,  -  Af„. 
I  t  Most  OBS’s  contain  low -density  flotation  for  retrieval 
in  their  upper  portion  and  more  dense  disposable  anchoring 
material  near  the  bottom-coupling  interface.  This  density 
gradient,  in  addition  to  producing  torques  from  equal  ampli¬ 
tude  horizontal  inputs,  will  modify  the  response  from  waves 
having  different  horizontal  amplitudes  in  the  water  and  in  the 
sediment  (for  example,  horizontally  propagating  comprcs- 
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sional  waves  or  horizontally  polarized  shear  wavesi  The  verti¬ 
cal  component  is  less  complicated,  since  signal  amplitude  is 
continuous  across  the  water -sediment  boundary  and  most 
OBS's  have  appreciable  symmetry  about  a  vertical  axis  How¬ 
ever.  as  pointed  out  by  tiurniany.  lateral  variations  within  the 
bottom  under  the  OBS  could  affect  this  component,  also 

The  response  of  an  OBS  to  a  vertical  seismic  input  was 
given  bv  Suitnn  ct  al  [19N14i]  and  /Wikoiir:  and  Prothcro 
|  1 9s  I  j  following  the  presentation  of  Sutton  et  al . 


where 

I  amplitude  of  OBS  in  inertial  space,  that  is.  input  to 
geophones ; 

amplitude  of  incident  wave. 

C  (Af,  ~  \f„  -  Af,t  Af,*: 

/i,  damping  coefficient  of  bottom  coupling. 
c\  resonant  frequency  of  bottom  coupling. 

s  transform  variable  Is  -  /<■»  for  sinusoidal  input I. 

Af,  instrument  mass: 

Af  „  mass  of  displaced  water: 

Af,  mass  of  displaced  sedimcm  I  usually  negligible): 

Af*  inertial  mass  of  instrument  plus  water  and  sediment 
that  move  with  it. 

Although  derived  for  vertical  motion,  when  interpreted  cor- 
recilv  the  equation  can  be  used  for  other  motion  components 
as  well.  The  transfer  function  contains  one  or  two  zeros,  de¬ 
pending  upon  (he  coupling  coefficient  (  two  zeros  if  ( '  *  l 
and  one  zero  if  C  -  I  The  case  (  -  I  implies  that  Af,  ~ 
\f,  -  0  and  A/,*  Af,:  that  is.  no  displaced  water  or  sediment 
and  no  additional  inertia  from  water  and  or  sediment  This 
form  of  the  equation  has  been  used  by  earlier  workers  to 
describe  the  motion  of  a  geophone  on  a  compliant  surface 
| e  g..  Solar.  1978],  The  zero  in  the  transfer  function  in  this  case 
is  related  to  the  radiauon  damping  from  the  differential 
mouon  between  the  geophone  and  the  soil  The  transfer  func¬ 
tion  is  also  consistent  with  the  relation  given  above  for  the 
moiton  of  a  sphere  in  a  fluid  In  the  absence  of  a  botiom- 
related  spring  constant. . ■»,  0.  and  when  ('  0  (that  is.  neu¬ 

tral  buoyancy.  Af,  -  Af.).  then  /  -  7.s  tl  ~  I  )  Additionally, 
when  < '  -  2.  then  /  -  A/,  il  I  =  As  mentioned  earlier, 

this  is  the  limiting  condition  when  \f,  -  Af,  -  0  and  Af,*  - 
W„  2.  \f„  2  is  the  appropriate  "virtual”  mass  to  be  added  to 
the  inertial  mass  of  a  sphere  moving  in  a  fluid  |  Batchelor. 
1967], 

(iumwim  [1984]  refers  to  Nutton  ef  a/  ( l9Klu.  b]  and  to 
/chkont:  and  Prothcro  [1981]  and  states  (hat  fin  lhc.se 
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papers)  data  from  the  Lopez  Island  OBS  experiment  ..  indi¬ 
cate  that  the  restoring  forces  due  to  buoyancy  are  2  or  more 
orders  of  magnitude  less  than  the  effective  spring  constant  of 
the  sediments.”  This  is  correct  and  means  that  the  modifi¬ 
cation  to  the  coupling  frequency  ojc  by  gravitational  restoring 
forces  on  the  OBS's  tested  at  Lopez  is  minor,  but  it  does  not 
mean  that  all  effects  of  buoyancy  are  negligible.  Table  V  of 
Sutton  et  al.  [1981b]  lists  values  of  Cm,x,  which  is  the  coupling 
coefficient  defined  above  with  M,  substituted  for  M*  in  the 
denominator  (that  is.  not  including  the  inertial  contribution  of 
the  water  and/or  sediment).  The  computed  values  of  CmiI  for 
operational  OBS’s  listed  in  Table  V  range  from  0.11  to  0.53; 
actual  values  for  C  are  still  smaller,  and  the  effect  of  buoyancy 
is  hardly  negligible.  For  neutral  buoyancy,  (7  =  0,  the  transfer 
function  predicts  perfect  coupling,  /  =  Z,,  independent  of  the 
coupling  frequency  and  damping. 

Garmany  [1984]  states  further 

The  fundamental  difference  between  the  present  work  and  pre¬ 
vious  work  on  OBS  coupling  is  that  the  author  believes  that  the 
motion  of  the  earlh  is  communicated  to  ihe  OBS  by  ihe  effective 
spring  constant  of  the  sediments  and  not  the  far  less  significant 
forces  of  buoyancy. 

He  suggests  that  proper  interpretation  of  the  position  vari¬ 
ables  used  by  Sutton  et  al.  [198  la]  would  indicate  that  buoy¬ 
ancy  had  negligible  effect  on  coupling  response. 

The  validity  of  the  transfer  function  given  above  is  difficult 
to  verify  with  field  experiments.  However,  Sutton  and  Duenne- 
bier  [1985]  reproduced  its  general  predictions  in  model  exper¬ 
iments  on  a  shaking  table  with  the  package  both  in  air  and 
submerged  in  water  and  for  both  vertical  and  horizontal 
inputs  and  outputs.  A  transfer  function  of  the  same  form  can 
also  be  used  for  a  rocking  mode  response  to  a  horizontal 
input  [Hsieh.  1962]. 

Garmany  [1984]  tested  his  procedure  using  OBS  data  and 
was  able  to  reduce  spectral  peaks  presumably  caused  by  cou¬ 
pling.  Unfortunately,  the  signals  resulting  from  this  procedure 
cannot  represent  “true  particle  motion"  as  they  are  obtained 
using  an  incorrect  parameterization.  Since  Garmany  assumes 
that  the  coupling  acts  as  an  all-pole  filter,  he  appropriately 
models  the  discrete  time  series  as  an  autoregression  (AR). 
However,  as  we  have  shown,  the  correct  response  contains 
both  poles  and  zeros,  so  the  correct  parameteric  time  series 
model  is  not  an  autoregression  but  an  autoregressive  moving 
average  (ARMA).  It  is  well  known  that  if  the  moving  average 
IMA)  part  of  an  ARMA  process  is  of  finite  order,  then  mod¬ 
eling  the  process  as  an  autoregression  will  require  an  infinite 
order  [Box  and  Jenkins .  1970,  chapter  3).  Hence  Garmany's 
treatment,  fitting  three-channel  AR  models  of  order  4  or  8 
only  (that  is.  12  or  24  poles),  can  never  be  expected  to  perform 
well.  A  dramatic  improvement  in  performance  could  be  ex¬ 


pected  if  an  AMR  A  model  were  fit  to  the  data  rather  than  a 
pure  AR  model.  If  our  theory  is  correct  (and  the  shake  table 
results  seem  to  confirm  it),  then  an  ARMA  (2.2)  (a  second- 
order  autoregressive,  sccond-ordcr  moving  average)  is  the  cor¬ 
rect  model  to  fit 

We  should  warn  that  any  parametric  model  fitting,  be  it 
AR,  MA.  or  ARMA,  is  designed  to  explain  correlations  in  the 
time  series.  If  such  a  parametric  model  is  used  as  a  basis  for 
filtering,  the  output  will  be  optimally  uncorrelated  (white)  for 
the  model  chosen.  Hence  such  procedures  have  the  effect  of 
reducing  all  spectral  peaks  and  will  reduce  those  related  to 
geology  as  well  as  those  caused  by  coupling.  Garmany  ac¬ 
knowledges  this  possibility. 

A  filter  obtained  by  fitting  a  parametric  time  series  model  to 
OBS  data  may  prove  valuable  in  reducing  the  effects  of  cou¬ 
pling,  but  it  cannot  be  a  substitute  for  better  OBS  design.  We 
agree  with  Garmany  that  the  value  of  proper  engineering  of 
OBS's  is  still  considerable  and  that  any  design  which  reduces 
the  in  situ  resonances  of  a  seismometer  will  be  a  significant 
advance.  It  is  clear  that  in  OBS  design  we  must  strive  for 
high-coupling  frequency,  near-critical  damping,  low  direct 
water  coupling,  and  low-coupling  coefficient  C  [Suffon  and 
Duennehier,  1985]. 
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